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ABSTRACT 


The  velocity  and  attenuation  of  longitudinal  sound 

4 

propagating  m  single  crystals  of  hep  He  at  a  molar 

3 

volume  of  17.4  cm  were  measured  between  0.25  K  and  melting, 
and  at  frequencies  of  5,  15  and  25  MHz.  All  samples  showed 
an  extremely  large  dispersion  and  in  most  cases  a  large 
velocity  anomaly  extending  from  between  1  and  1.8  K  down 
to  the  lowest  temperatures.  At  higher  frequencies  the 
anomaly  was  reduced  in  magnitude  and  it  shifted  to  lower 
temperatures.  The  anomaly  was  usually  positive,  but  in 
two  cases  it  was  negative  and  in  two  other  cases  was  not 
present  at  all.  The  attenuation  was  0.01  cm  ^  or  less 
near  melting  and  increased  up  to  100  times  by  the  lowest 
temperature,  except  for  the  two  crystals  without  velocity 
anomalies;  these  showed  no  dramatic  rise  in  attenuation. 

The  attenuation  increased  with  frequency.  Several  theories 
involving  second  sound,  zero  sound,  dislocations  and 
superfluidity  are  invoked  in  order  to  explain  the  effects. 
None  was  successful. 
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CHAPTER  1 


INTRODUCTION 


Solid  helium  was  originally  studied  with  the 

expectation  that  it  would  prove  to  be  the  simplest  and 

most  ideal  of  all  solids.  For  some  physical  properties 

these  hopes  have  been  at  least  qualitatively  realized, 

but  in  other  respects  the  behaviour  of  solid  helium 

deviates  dramatically  from  that  of  more  normal  solids. 

Most  of  these  unusual  effects  can  be  attributed 

to  the  quantum  nature  of  helium.  Eecause  the  atoms  are 

very  light,  and  interact  only  through  a  weak  Van  der  Waals 

potential,  they  possess  a  very  large  zero  point  energy  and 

experience  such  large  excursions  from  their  lattice  sites 

that  the  wave  functions  of  neighbouring  atoms  overlap  to 

a  significant  extent.  The  symmetry  of  the  wave  functions 

3 

then  becomes  important  so  that  liquid  Pie  ,  a  fermi  particle 

(2  protons,  1  neutron  and  2  electrons)  has  markedly  dif- 

4 

ferent  properties  from  liquid  He  ,  a  bose  particle 
(2  protons,  2  neutrons,  2  electrons).  The  localization 
of  the  atoms  on  lattice  sites  in  the  solid  reduces  the 
overlap  so  that  the  differences  between  the  two  isotopes 
are  more  subtle  in  the  solid  state. 

Since  it  has  a  lower  boiling  point  than  any  other 
substance,  helium  can  be  purified  to  an  extremely  high 
degree  with  simple  techniques.  Only  disorder-type 
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. 
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defects  will  then  be  present.  As  well,  single  crystals 
of  solid  helium  are  relatively  easy  to  grow,  and  it  is  a 
dielectric.  For  these  several  reasons,  research  on  solid 
helium  is  largely  free  of  the  many  spurious  effects  which 
plague  the  study  of  other  more  complex  solids.  Transport 
properties  in  particular  may  be  observed  under  the  nearly 
ideal  conditions  which  are  difficult  and  expensive  to 
attain  in  other  crystals.  A  lightly  damped  second  sound 
mode  may  be  propagated  while  the  thermal  conductivity 
exhibits  the  pronounced  umklapp  maximum  and  well-defined 
Poiseuille  flow  region  typical  of  very  pure  dielectrics. 

Solid  helium  has  an  extraordinarily  high  compres¬ 
sibility,  enabling  observations  to  be  taken  over  a  wide 
range  of  molar  volumes  (the  volume  is  halved  with  the 
application  of  1500  bar  of  pressure) .  Of  thermodynamic 
interest  is  the  fact  that  solid  helium's  high  compressibi¬ 
lity  allows  its  properties  to  be  measured  at  an  essentially 
constant  volume. 

The  present  work  may  be  considered  in  part  to  be 
the  third  stage  of  a  long  term  project  to  measure  the 
temperature,  volume  and  frequency  dependence  of  the  attenu¬ 
ation  of  sound  in  solid  helium.  I  believe  the  data  con¬ 
tained  in  this  work  to  be  the  first  reliable  measurements 
of  attenuation,  particularly  at  higher  temperatures. 

The  most  important  goal  of  this  work  is  to  examine 
in  greater  detail  the  anomalous  behaviour  of  the  sound 
velocity  which  was  discovered  in  this  laboratory  by  Franck 
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and  Hewko  (1973)  and  further  investigated  by  Wanner  and 
Mueller  (1974)  and  Wanner  et  al.  (1976).  The  original  aims 
were  to  observe  the  anomaly  with  much  greater  precision 
and  to  extend  the  measurements  to  lower  temperatures.  The 
high  quality  of  the  ultrasonic  transducers  used  enabled 
me  to  successfully  excite  and  observe  the  third  and  fifth 
harmonics  of  the  fundamental  frequency.  This  allowed  not 
only  the  dispersion  in  solid  helium  to  be  measured  for  the 
first  time  in  the  long  wavelength  limit,  but  also  the 
frequency  dependence  of  the  attenuation  and  the  velocity 
anomaly . 

Preliminary  results  of  this  work  have  already 
been  published  elsewhere  (Calder  and  Franck  (1977) ) . 


CHAPTER  2 


THEORY 

For  the  study  of  classical  lattice  dynamics  and 
sound  propagation  I  have  found  particularly  useful  books 
by  Maradudin  et  al.  (1971)  ,  Born  and  Huang  (1954) , 

Musgrave  (1970)  and  Reissland  (1973),  with  supporting 
material  from  Nye  (1957) .  There  are  available  numerous 
useful  reviews  of  quantum  lattice  dynamics  including  those 
by  Guyer  (1969) ,  Koehler  (1975)  and  Glyde  (1976) . 

The  most  recent  reviews  of  the  properties  of  solid 
helium  are  the  works  of  Wilks  (1967)  ,  Keller  (1969)  , 
Trickey  et  al.  (1972)  and  Glyde  (1976) .  Theories  of  the 
various  contributions  to  the  temperature  dependence  of  the 
sound  velocity  and  attenuation  can  be  found  in  Bhatia 
(1967),  Truell  et  al.  (1969)  and  especially  Maris  (1971). 

2 . 1  Classical  Lattice  Dynamics 

It  is  obvious  from  numerous  experiments  like  the 
present  one  that  sound  does  propagate  in  bulk  solid  helium 
in  the  long  wavelength  limit.  With  this  observation  in 
mind,  we  can  confidently  go  ahead  and  use  the  relations  of 
continuum  mechanics  as  developed  in  section  2.3  in  order 
to  relate  the  velocity  of  sound  to  such  macroscopic 
properties  of  the  crystal  as  the  elastic  constants.  How¬ 
ever  in  order  to  use  acoustic  modes  effectively  as  probes 
of  the  properties  of  a  solid  on  the  atomic  scale,  we  need 
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a  microscopic  theory  of  lattice  dynamics.  Although  clas¬ 
sical  lattice  dynamics  are  totally  inadequate  to  describe 
solid  helium,  even  with  anharmonic  corrections,  a  brief 
outline  still  serves  as  a  useful  introduction  to  the 
methods  required  when  dealing  with  solid  helium. 

Each  atom  sits  in  the  potential  well  created  by 
all  other  atoms  in  the  crystal.  If  the  atoms  make  only 
very  small  excursions  from  their  equilibrium  positions, 
then  the  potential  may  be  expanded  into  the  harmonic  form 


V 


V  + 
o 


2 

9  V 


9u 


s£ 


8u 


s  '  £  1 


0 


u 


j’ 

s  '  £  ' 


(2.1-1) 


where  there  is  summation  over  repeated  indices  and  u^ 

represents  the  scalar  displacement  of  atom  s  in  unit 

cell  £  in  the  direction  j .  Higher  order  terms  may  be 

added  to  (2.1-1)  to  give  anharmonic  corrections.  The 

total  kinetic  energy,  for  atoms  of  mass  M  is: 

s 


T 


(2.1-2) 


The  equations  of  motion  can  then  be  easily  found  from 
Lagrange's  equations  giving, 


' 
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where 


D 


D  J 

s  £  ;  s  '  £ 


82V 


Su-1  „  8u^  ,  , 

s£  s  £ 


(2.1-4) 


0 


and  henceforward  there  is  no  summation  over  s  only. 

Assuming  only  central  forces  are  present,  and 
using  the  Bloch  condition, 


we  find, 


Hs£(t)  =  el-‘"  Uso(t)' 


M  ii  (t)  =  -  D  ,  (q)  *u  ,  (t)  , 

s-sq  =ss'  _  -s'q 


(2.1-5) 


(2.1-6) 


where  q  is  the  wave  vector  and. 


D  , <9)  =  l  E  . .  -lq,h 


=ss  ’  ’-±'  t  =s s  1  ; h 
h 


e  -  - 


(2.1-7) 


Assume  a  plane  wave  solution: 


u  (t)  =  u  e 

-sq  -sq 


imt 


(2 . 1-8) 


so  we  have  finally, 


[D  ,(q)~  M  co  (q)  6  ,I]*u  ,  =  0 

=ss'  2  s  -  ss'=J  -s'q 


(2.1-9) 


where  6  ,  is  the  Kronecker  delta  and  I  the  identity 

s  s 

matrix.  One  can  now  find  co  as  a  function  of  q  by 
solving  the  secular  equation, 
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l£ss,(q)  "  Msw2(3),5ss'il  =  °-  (2.1-10) 

This  equation  has  3n  solutions  where  n  is  the  number  of 
particles  in  one  unit  cell  and  the  three  refers  to 
j  =  1,  2,  3,  the  three  polarizations  of  the  vibrations. 

The  group  velocity  is  defined  as, 


c  . 
-J 


V  •  m  .  (q)  , 

-q  J  - 


(2.1-11) 


but  usually  what  one  measures  is  the  phase  velocity  , 


defined  by. 


°°j  (5)  =  -j  (q)  'q  * 


(2.1-12) 


In  the  case  where  there  is  no  dispersion,  the  phase  velo¬ 
city  is  identical  with  the  group  velocity. 

In  the  long  wavelength  limit  (i.e.  for  small  q) , 

D  , (q)  in  (2.1-10)  can  be  expanded  in  powers  of  q 

S  S  —  “■ 

(Maradudin  et  al  (1971))  and  (2.1-10)  can,  after  some 
lengthy  algebra,  be  written  in  the  form, 


(C 


a$  Ay 


+  C'  , 0  )  q,  q 
aA3y 


-  p  co  .  6  0 

]  aB 


=  o. 


(2.1-13) 


where  C  and  C'  are  complicated  functions  of  the 

derivatives  of  D  , (q) .  Equation  (2.1-13)  provides  a 

s  s 

link  between  lattice  dynamics  and  the  continuum  model 
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discussed  in  section  2.3. 

We  can  gain  further  insight  into  classical  lat¬ 
tice  dynamics  by  forming  the  Hamiltonian  of  the  system  in 
the  harmonic  approximation  from  (2.1-1)  and  (2.1-2)  and 
diagonalizing  it  to  give: 

H  =  h  £  [  |  p  (k  j  )  |  2  +  m2  (k  j  )  |q  (kj  )  |  2]  (2.1-14) 

kj 

where 

u“£  =  £'(g£a;kj)  q(kj)#  (2.1-15) 

u“^  =  B,  ( s  £  a ;  k  j  )  p(kj),  (2.1-16) 

2  rv  ry  * 

and  03  (kj)  =  B,*  (s£a;  kj  )  £  (s  '  £  '  a  ;  k '  j  ' )  6R_k  ,  , 

(2.1-17) 

with  £(s£a;kj)  =  M  ^  e  (kj)e  -s  .  (2.1-18) 

"—P  s  ct 

Using  (2.1-18)  in  (2.1-17)  along  with  the  Bloch  theorem 
once  again  we  find: 

[D  (k  j  )  -  mo)2  (k  j  )  I]  u  (k  j  )  =  0,  (2.1-19) 

as  in  equation  (2.1-9). 

Equation  (2.1-14)  then  gives  the  Hamiltonian  for 
a  harmonic  oscillator  so  that  the  quantized  lattice  modes, 
or  phonons,  obey  Bose-Einstein  statistics. 
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2 . 2  The  Lattice  Dynamics  of  Solid  Helium 

In  most  normal  crystals  the  harmonic  approximation 
used  in  equation  (2.1-1)  is  a  good  one.  As  long  as  the 
kinetic  energy  of  a  particle  is  only  a  small  fraction  of 
the  depth  of  the  potential  well,  then  the  particle  remains 
very  near  to  the  bottom  of  the  potential  well,  where  (2.1-1) 
is  valid.  Small  excursions  from  (8V/3u)  =  0  can  usually 
be  accounted  for  with  small  anharmonic  corrections  to 
equation  (2.1-1). 

The  interatomic  forces  for  helium  are  the  Van  der 
Waals  forces,  which  are  very  weak  for  the  tightly  bound 
electronic  shell  around  helium  atoms.  A  simple  approxi¬ 
mation  is  the  Lennard-Jones  potential: 

V (r)  =  4e  [o/r)12  -  (o/r)6],  (2.2-1) 

as  shown  by  the  dashed  curve  of  figure  1.  From  the 
virial  coefficients  of  helium  gas,  it  is  found  that  the 
helium  potential  has  a  very  shallow  well  of  e  =  10.2  K 
and  a  very  extended  hard  core  with  a  =  2.556A.  If 
(2.2-1)  is  used  in  equation  (2.1-10),  then  the  sound 
velocity  predicted  by  (2.1-10)  is  four  times  larger  than 
that  found  experimentally.  Classical  lattice  dynamics 
also  predict  a  molar  volume  too  small  by  half  and  a  com¬ 
pressibility  thirty  times  too  small.  On  the  other  hand, 
heat  capacity  and  thermal  conductivity  behave  quite 
normally,  so  that  some  sort  of  long  range  correlations 


V ( r )  (K) 


10 


r  (A) 


Figure  1  Interatomic  potential  well  in  solid  helium-4  at  17.4  cm^/mole. 
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must  be  present. 

At  P  =  0,  T  =  0,  the  kinetic  energy  of  a  solid 

arises  solely  from  the  zero  point  energy  e.  .  In  the  case 

of  helium  e  -  V(r  ),  where  r  =  2 ^^•o,  so  that  solidi- 
o  o  o 

fication  will  only  take  place  with  the  application  of  a 
substantial  pressure  (at  least  25  bar) .  If  the  helium 
atoms  were  to  sit  at  the  bottom  of  the  potential  well 
created  by  neighbouring  atoms,  then  the  large  excursions 
resulting  from  the  zero  point  dynamic  energy  would  bring 
neighbouring  atoms'  hard  cores  into  close  contact.  Thus 
the  atoms  move  farther  apart  so  that  we  have  the  situation 
shown  in  figure  1.  In  this  model,  the  atoms  sit  at  local 
potential  maxima  so  that  their  vibrational  frequencies 
0)  v  27r(V"(r)/m)  2  become  imaginary,  as  first  calculated  by 
DeWette  and  Nijboer  (1965) .  Thus  we  have  a  complete 
breakdown  of  classical  lattice  dynamics. 

Any  successful  lattice  dynamical  theory  must  take 
into  account  not  only  the  long  range  correlations,  or 
phonons,  but  also  the  short  range  correlations  of  neigh¬ 
bouring  atoms  which  move  in  a  cooperative  manner  in  order 
to  minimize  wave  function  overlap. 

The  method  which  allows  the  simplest  physical 
interpretation  (Nosanow  (1964,  1966),  Fredkin  and  Wert- 
hamer  (1965) ) ,  assumes  the  wavefunction  proposed  by 


Jastrow  (1955) : 
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¥  =  n<#> .  (r  . -£  .  )  IT  f  (r.-r  .)  , 
.11  1.  .  vi  i 

i  i<  j  J 


(2.2-2) 


where  <j>^(r^)  are  the  single  particle  wave  functions  local¬ 
ized  at  lattice  sites  while  f(r.-r.)  is  a  Jastrow 

i  i  j 

function  describing  the  short  range  correlation  between 
particles  i  and  j.  A  useful  analytic  form  for  f (r)  is: 


f (r)  =  e 


-KV (r) 


(2.2-3) 


where  K  is  a  variational  parameter.  Now  (p  is  so 
close  to  a  Gaussian  that  we  let 

-Ar . 2/2 

(p.  (r.)  =  e  1  (2.2-4) 

Yi  l 

and  if  ¥  is  used  in  the  first  term  of  a  cluster  expan¬ 
sion  for  the  energy  E  (i.e.  in  a  variational  approach). 


Hi 

2 


(2.2-5) 


then  we  find 


E  =  kh 


£  00,  ,  +  £ 

i  A  k ,  A  .  . 

k, A  -  i< j 


f$2  V 


eff 


(rij) 


(2.2-6) 


where 


$  =  II  <J>.  (r .  )  • 

.  Yi  i 
i 


(2.2-7) 


Veff  (r__)  is  an  effective  potential  given,  for  two-body 
interactions,  by. 


Veff(rij) 


f  (rij)[V(rij> 


V2V  (r .  .)  ]  . 
2m  13 


(2.2-8) 


Then  K  and  A  are  varied  to  minimize  E. 

Now  the  long  range  correlations  must  be  properly 
accounted  for  by  using  self-consistent  phonon  theory 
(Nosanow  and  Werthamer  (1965)  e  De  Wette  et  al.  (1967)). 
This  is  essentially  a  Hartree  method  in  which  the  fre¬ 
quency  00  becomes  a  variational  parameter  so  we  set 


dE _  -h  „  ,  .  d$2 

9a3kX  4  i<j  eff  r±3  9t\x 


0 


(2.2-9) 


Then  00^  are  given  by  the  determinant  of  the  dynamical 
matrix 


(k) 


(e 


-ik  • 


i .  e . 


0. 


(2.2-10) 

(2.2-11) 


Equation  (2.2-11)  resembles  equation  (2.1-10)  very 
closely.  Only  the  dynamical  matrix  is  defined  differently 
than  in  equation  (2.1-7)  so  that  the  second  derivatives 
of  the  potential  are  now  averaged  in  order  to  account  for 
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long  range  correlations,  while  an  effective  potential  is 
used  to  allow  short  range  correlations.  With  these  re¬ 
normalized  force  constants  we  have  a  quasiharmonic  result 
in  which  we  can  still  speak  of  phonons  as  good  excitations 
in  solid  helium. 

Several  improvements  have  been  made  by  using 
higher  order  terms  in  the  cluster  expansion  of  equation 
(2.2-5)  (Hetherington  et  al.  (1967)),  by  using  a  fully 
correlated  wave  function  rather  than  a  Gaussian  approxi¬ 
mation  (Koehler  (1967) )  and  by  including  higher  order 
anharmonic  terms  in  equation  (2.1-1)  (Koehler  (1968), 
Koehler  and  Werthamer  (1972) ) .  Other  recent  developments 
in  the  dynamics  of  quantum  crystals  are  discussed  by 

Glyde  (1976)  and  Koehler  (1975) .  Most  of  these  have  been 

3 

applied  only  to  He  . 


2 . 3  Wave  Propagation  in  Ideal  Elastic  Solids 

An  excellent  reference  work  upon  which  I  have 
relied  heavily  for  the  contents  of  this  section  is  the 
book  by  Musgrave  (1970). 


If  u^(Xj)  is  the  displacement  in  the  i  direction 
as  a  function  of  x^ ,  the  position  within  a  continuous 
elastic  medium,  then  for  small  values  of  the  displacement 


we  define  the  strain  tensor  as 


+ 


(2.3-1) 


£■£ 
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As  well,  the  stress  a.  .  is  defined  as  the  force  per  unit 

ID 

area  acting  in  the  i  direction  upon  a  surface  with 
normal  in  the  j  direction.  For  small  strains  Hooke's 
law  holds  so  that  in  general: 


Qij  Cijk£  ek£  ' 


(2.3-2) 


where  the  elements  of  the  tensor  c.  „  are  called  the 

ijk£ 

elastic  stiffness  constants,  and  from  dynamic  considera¬ 
tions  , 


ci  jk£ 


32V 


8  £  .  .  8  £,  0 

13  k£ 


(2.3-3) 


By  symmetry  arguments  the  number  of  independent  elastic 
constants  can  be  substantially  reduced.  Using  the  Voigt 
notation  where  c^j^  becomes  a  6x6  matrix  and  index  pairs 
become  single  indices  according  to  the  prescription 


rH 

+ 

1 — 1 

1 — 1 

22  ->  2 

,  33  -* 

3,  23 

or  32 

-*  4, 

31  or  13 

■>  5 ,  12 

or  21  -> 

6  ,  c  .  . 
ID 

for  a 

crystal  with 

hexagonal  symmetry  is : 

r 

C11 

c12 

c13 

0 

0 

* 

0 

C12 

C11 

C13 

0 

0 

0 

c  .  .  = 

C13 

c13 

C33 

0 

0 

0 

(2.3-4) 
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0 

0 

0 

c.. 

44 

0 

0 

0 

0 

0 

0 

C44 

0 

0 

0 

0 

0 

0  h 

(C11_C12) 
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The  force  on  an  infinitesimal  volume  is 


3o  .  . 

=  m 


3x  . 
1 


pui , 


(2.3-5) 


so  that  upon  substituting  (2.3-1)  and  (2.3-2)  into 
(2 . 3-5) ,  we  find, 


32un 


c .  . 


i jk£  3xn  3x  .  °Ui  ^  * 
i  j 


(2.3-6) 


If  we  assume  a  plane  wave  form  for  u. 


uk  =  Ak  exp  i  (q-x  -  out)  , 


(2.3-7) 


then 


cijk£  qjqil  -  pw  6ikl=  0 


(2.3-8) 


is  the  resulting  secular  equation  from  which  w(q)  may  be 
obtained.  Equation  (2.3-8)  is  identical  to  equation 
(2.1-10).  Several  authors  including  Zener  (1936),  Gold 
(1950)  and  Musgrave  (1954)  have  solved  equation  (2.3-8) 
explicitly  for  the  hexagonal  case.  They  find  that  the 
velocity  surfaces  have  cylindrical  symmetry  about  the 
c-axis  so  that  all  velocities  may  be  expressed  in  terms  of 
a  single  angle  y  between  the  c-axis  and  the  direction  of 
propagation.  There  are  three  mutually  perpendicular 
polarizations:  a  pure  transverse  mode  T^,  a  quasi- 
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traverse  (i.e.  predominantly  transverse)  mode  T ^  and  a 
quasi-longitudinal  mode  L.  These  phase  velocities  v  are 
given  by 


pvr 


pvr 


PVT 


2 

2 


where 


=  ^ (cll“cl2) sin2y  +  c44cos2y  (2.3-9) 

2  2 

=  hi  (c11+c44)  sin  y  +  (c33+c44)cos  y  —  cf>  (y )  ] 

(2.3-10) 

2  2 

=  (c11+c44)  sin  +  (c33+c44)cos  y  +  <f>(y)], 

(2.3-11) 

,  2  ,  .  ,  N  2  .  4  ,  ,  ,2  4 

4>  (Y)  =  ( cn~c4 4 )  Sln  Y  +  ^c33~c44  ^  cos  Y  + 


+  2sin2Ycos2y  [(ci:L-c44)  (c44-c33)  +  2  (c13+c44 )  2  ]  . 

(2.3-12) 

3 

The  velocities  for  a  molar  volume  of  17.4  cm  are  illus¬ 
trated  in  figure  2,  using  the  elastic  constants  calculated 
in  appendix  1. 

Without  loss  of  generality  we  may  set  the  azimu¬ 
thal  angle  equal  to  90° ,  so  that  the  polarizations  are 
given  by 


=  d,  0,  0) 

il 


(2.3-13) 


0 


(2.3-14) 


006 
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Figure  2  Orientation  dependence  of  the  sound  velocity  at  17.4  cm3/7mole. 
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(pL,T  =  f  (PV2,cij  ,y)/ [1  +  f2  (pv2,ci;.  ,y)  ]  ^ 


(pL,T  >3  =  I1_(PL,T  >22]!i  =  [1  +  f2<PV2,Cij,Y) 


(2.3-15) 


(2.3-16) 


2  /  Cll+C12 

where  f  (pv  ,  c.  .  ,y)  =  Jg  - — — —  )  tany  x 

13  \  c13  44 


x 


2  2  2  2 
pv  “c44sin  Y  “  tc33-2  (ci3+c44)  /  (cH+cl2 )  ]  cos  Y 


2  2  ,  ,  v  •  2 

pv  “c44cos  Y  “  *s(c11-c12)sin  y 


(2.3-17) 


The  angle  6  which  the  direction  of  polarization  makes 
with  the  direction  of  pure  transverse  or  longitudinal 
vibration  is  shown  in  figure  3  as  a  function  of  orienta¬ 
tion  y  for  the  and  L  modes  at  a  molar  volume  of 
17.4  cm  .  The  angle  between  the  polarization  vector  and 
the  c-axis  is  y-6,  since  the  c-axis,  the  polarization 
vector  and  the  propagation  direction  are  coplanar. 

In  general  the  direction  of  propagation  of  sound 
is  not  perpendicular  to  the  transducer  but  is  at  some 
angle  A  from  the  normal  direction: 


cos  A, 


1  + 


r/  ,2  .2,2.2  2 

[(C11_C12)  "  4c44-*  Sin  YCOS  Y 


r  ,  v  2  .  2  ,  . 

[ (cii-ci2^  Sin  Y  +  4c 


44 


2  2.2 

cos  ] 


-h 


(2.3-18) 


Deviation  of  Beam  Direction 
and  Polarization  8  at  17.4 
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Figure  3  Orientation  dependence  of  the  beam  deviation  (A)  and  polarization  (§)  at  17.4  cm3/mole. 
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(2.3-19) 


where  are  the  polarization  directions.  The  beam 
deviations  for  a  molar  volume  of  17.4  cm  are  given  in 
figure  3. 


2 . 4  Thermodynamic  Quantities 

Several  thermodynamic  quantities  of  interest  are 
closely  related  to  the  sound  velocity. 


Compressibility  and  Elastic  Constants 

The  compressibility  can  be  defined  at  constant 
temperature  (isothermal)  or  constant  entropy  (adiabatic) 
as  (Nye  ( 19  57  )  )  , 


T  ,  T  .  T  ,  0  T 
C11  C12  4c13  2c33 

.  T  T  ,  T  0  ,  T  .  2 .  ' 

(C11+C12)C33  "  2  C1 3 


s  s  .  s  .  n  s 

cll  C12  4c13  2c33 
-S^S.S  0  ,  S  x  2  ' 

Cll+C12  C33  2  C13 


(2.4-1) 


(2.4-2) 


T  s 

where  c.  .  and  c.  .  are  elastic  constants  defined 
iD  ID 

isothermally  and  adiabatically .  Using  some  simple  thermo¬ 
dynamic  relations,  the  difference  between  adiabatic  and 
isothermal  elastic  constants  is  easily  shown  to  be 


22 


(o-j-c-.)/c-j  *  (Cp-Cv)/Cv  <1.  T4,  (2.3-3) 

where  Cp  and  are  the  heat  capacities  at  constant 

4 

pressure  and  constant  volume.  For  the  case  of  He  ,  where 
the  ratio  of  the  lattice  parameters  c/a  is  independent 
of  volume,  then  as  shown  by  Franck  and  Wanner  (1970), 


CU  +  c12  =  c33  +  c13 


(2.4-4) 


so  the  compressibility  becomes 


k  3/  (c33  +  2c13) • 


(2.4-5) 


Debye  Temperature 

The  specific  heat  of  a  Debye  solid  at  low 
temperatures  is  (Maradudin  et  al.  (1971) ) , 


0  is  a  function  of  volume  only  for  the  ideal  Debye  solid, 
but  for  real  solids  it  is  also  a  slowly  varying  function  of 
temperature.  The  Debye  temperature  at  T  =  0  is  related 
to  the  sound  velocities  (not  necessarily  at  T  =  0)  by 


(2.4-7) 
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where  the  Debye  velocity  is: 


r  r  ,  -3,  -3,  -3  >  dfi  , 

vD  =  t/(vL  +  vTi  +  vT^  )  127  ] 


Griineisen  Constant 


-1/3 


(2.4-8) 


The  Griineisen  constants  are  in  general  defined 


by 


Y 


kj 


3£nco,  . 

]S  3 

3  £nV  * 


(2.4-9) 


Note  that  since  the  dynamical  matrix  in  (2.1-10)  is 

independent  of  volume,  y^  will  only  be  non-zero  for  an 

anharmonic  solid.  Now  from  (2.4-7),  (2.4-8)  and  (2.1-12), 

if  y,.  is  independent  of  k  and  j,  then  we  have  only  one 
_3  ~ 

Griineisen  constant  given  by 


3  £n0 


Yo  3  £nV 


o 


(2.4-10) 


If  a  system  can  be  described  by  a  reduced  equation  of 
state , 

(MV)  =  <J)  (V)  -f  [T-(f>  (V)  ]  ,  (2.4-11) 

where  usually  4>  (V)  =  0  ^ ,  then 


3£n6  V  3P 

3£nV  C  3T 

v 


(2.4-12) 


v 


■ 
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Second  Sound 

The  propagation  of  second  sound  may  be  understood 
on  an  elementary  level  by  using  the  method  of  Ward  and 
Wilks  (1952)  with  the  Boltzmann  equation  of  Callaway 
(1959)  : 


-  p-Vf 
P  -  - 


(2.4-13) 


t  and  tkt  are  the  total  relaxation  times  for  resistive 
and  momentum-conserving  phonon  processes.  The  resistive 
processes  relax  the  distribution  to  fQ  and  the  conserva¬ 
tive  forces  relax  it  to  f ^ .  Defining  the  energy  density 
as : 

E  =  /cpfd^p,  (2.4-14) 


it  is  easily  shown  that. 


afl  +  1_  M  -  si  v2 

_  2  t  9t  3 
9t  u 


(2.4-15) 


Thus  second  sound  is  a  damped  oscillation  of  the  energy 
density.  We  assume  a  plane  wave  solution. 


E 


E  +  E  e1 
o  1 


(2.4-16) 


Then  the  real  and  imaginary  parts  of  the  wave  number 
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k  are : 


and 


(2.4-17) 


(2.4-18) 


If  (jot^&I,  then  k_^/k^l  and  we  have  diffusive  propagation. 

But  if  (jl) t  >>1,  then, 
u 


C2  =  w/kr  =  c 


(2.4-19) 


and 


/3  03  1 

Ot  ^  .  .  . 

2  2  C  GJT 

u 


(2.4-20) 


The  second  sound  velocity  has  been  calculated  for  several 
polarizations  of  phonons  in  an  isotropic  material  by 
Sussmann  and  Thellung  (1963)  to  be: 


c2  =  |  (I  vT3)/(£  v..5)  .  (2.4-21) 


c 2  has  also  been  calculated  for  the  more  general  aniso¬ 
tropic  case  by  Beck  (1975),  Banerjee  and  Pao  (1974), 

Hewko  (1973)  and  Hardy  (1970) . 

The  statistical  average  implied  by  the  relaxation 
time  approximation  of  equation  (2.4-13)  must  be  taken 
over  a  time  much  longer  than  the  total  relaxation  time. 


. 
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i.e.  ooU^1  +  tj(j1  )  1  <<  1. 


(2.4-22) 


Now  (wTy)  ^<<1,  so  (2.4-22)  is  satisfied  only  if 


COT  <<1 . 
N 


(2.4-23) 


Thus  we  have  a  window  for  second  sound  propagation,  as 
originally  shown  by  Guyer  and  Krumhansl  (1964)  : 


(2.4-24) 


ojtn<<1  <  <  go  t 


u 


2 . 5  Adiabatic  Sound  Propagation 

In  order  for  sound  to  propagate  adiabatically 
(in  quasi-equilibrium) ,  the  heat  which  flows  across  one 
wavelength  must  be  much  less  than  that  required  to  equal¬ 
ize  the  temperature  (Bhatia  (1967)). 


i.e.  Q1  <<  Q2 


(2.5-1) 


where  Q 


(2.5-2) 


(2.5-3) 


Therefore 


(2.5-4) 


Equation  (2.5-4)  means  that  sound  propagates  adiabati- 
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cally  at  high  temperatures  or  low  frequencies. 

From  a  microscopic  viewpoint,  a  sound  wave  in  an 
ideal  crystal  will  be  attenuated  and  shifted  in  velocity 
for  three  separate  reasons,  each  operative  under  a  differ¬ 
ent  set  of  conditions.  At  low  temperatures  the  sound 
waves  behave  like  a  beam  of  phonons  which  is  scattered 
by  the  thermal  phonons.  This  mechanism  will  be  discussed 
in  section  2.7.  At  higher  temperatures,  cot  <1,  a  longi¬ 
tudinal  wave  will  create  a  periodic  density  variation  so 
that  regions  of  compression  are  warmer  than  those  of 
rarefaction.  Heat  is  then  conducted  irreversibly  between 
these  regions  either  through  normal  diffusive  heat  con¬ 
duction  or  through  second  sound  propagation,  as  discussed 
in  section  2.6. 

Also  at  higher  temperatures,  when  oot<<1  (adia¬ 
batic  region) ,  the  phonon  gas  is  strained  by  the  passage 
of  the  sound  wave.  The  phonon  frequencies  will  be  locally 
modified  through  equation  (2.4-9)  and  then  will  relax 
towards  their  equilibrium  valve,  so  that  energy  is  dis¬ 
sipated.  Several  authors  have  investigated  this  mechan¬ 
ism,  including  Akhieser  (1939)  ,  Pomeranchuk  (1941) , 

Bommel  and  Dransfeld  (1960),  Woodruff  and  Ehrenreich  (1961) 
and  Maris  (1967,  1971) . 

In  this  temperature  range  a  Boltzmann  transport 
equation  approach  is  valid.  Assuming  only  umklapp  pro¬ 
cesses  are  important,  and  that  y,  .  is  a  constant  for  all 
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modes.  Woodruff  and  Ehrenreich  (1961)  find  the  attenuation 
is  given  by, 

Y^C  T 
7T  '  V 

a.  -  - -x —  .  cot  .  (2.5  —  5) 

J  pVv 


Maris  (1971)  has  calculated  the  adiabatic  velocity  v^  in 
this  region  to  be: 


v  -v 
a 


o 


_ 1_ 

2  pVv 


I 

kj 


Ys  (k  j  ) ho3  (kj  )  [n  (kj  ) 


+ 


x  (<Y^(kj)>  “  <Yc(kj)>2)> 

o  S 


C  T 
V 


x 


(2.5-6) 


where  n(kj)  is  the  phonon  occupation  number, 


Ys (kj) 
Ys (kj) 


Y  (kj)  - 


9  Y  (k j ) 

9^—  ‘  eY^J)*K6  ' 


1 


y6 

3co  (k  j  ) 


oo  (kj )  3e 


a8 


.  e  (KJ)  -K, 

a  - 


(2.5-7) 


(2.5-8) 


and  <f(kj)>  is  f(kj)  averaged  over  the  heat  capacities  of 

all  phonon  modes.  e(KJ)  is  the  polarization  of  the  sound 

/\ 

wave  KJ;  K  is  a  unit  vector  in  the  K  direction,  and  e  0 

-  -  -  ap 

is  the  strain  tensor.  Equation  (2.5-8)  is  similar  to 
(2.4-9).  For  the  case  of  an  isotropic  material  obeying 
a  reduced  equation  of  state,  (2.5-6)  reduces  to: 


. 


u. 


(2.5-9) 


v  -v 
a 


o 


_ L 

2pVv 


d£ny 

dJinV 


Including  one  higher  order  term  in  the  internal  energy  U 
(Maradudin  et  al.  (1971) ) ,  (2.5-9)  implies: 

Av/v  ^  aT4  +  bT6.  (2.5-10) 

2 • 6  Sound  Propagation  at  Intermediate  Temperatures 

At  this  point  the  various  phonon  scattering 
processes  must  be  examined  more  closely  (see  Reissland 
(1973)).  Resistive  processes  are  of  several  types.  Scat¬ 
tering  from  impurities  is  usually  insignificant  in  helium 
because  of  the  high  purity  obtainable.  A  discussion  of 
scattering  from  other  lattice  imperfections  such  as  dis¬ 
locations  and  vacancies  will  be  left  to  sections  2.8  and 
2.9.  Phonons  in  an  ideal  crystal  still  experience 
umklapp  scattering: 


(2.6-1) 


where  k^  and  k^  are  initial  and  final  phonon  momenta 
while  K  is  a  reciprocal  lattice  vector.  Energy  is  always 


conserved : 


1 ;y 


.  4 


4  J 
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For  an  n-phonon  process,  at  least  one  phonon  must  have 
I  IS  I  >  |K/(n-l)  |,  so  that  at  low  temperatures  umklapp 
processes  become  unimportant. 

If  K  =  0  in  (2.6-1),  then  we  have  normal  processes. 
For  three-phonon  processes  when  the  frequency  of  the 
sound  wave  03  <<  kT/fi,  then  (2.6-2)  becomes 


a) 1  =  co  +  oo  =  m+  (V.  <o)  •  Ak  . 

s  k  — 


(2.6-3) 


Using  the  group  velocity  v^ 
defined  by  (2.1-11)  and  (2. 
maximum  value  of 


and  the  phase  velocity 
1-12),  03 '  -  03  will  have 


a 


A 03  =  v  (k  j  )  I  •  I  Ak  I  , 

max  1 -g  - J  1  1  —  1 


(2.6-4) 


and  A  03  =  o)  =  v  •  I  k  I  . 

s  1 -p 1  1 -s 1 


(2.6-5) 


Therefore  |v  Ilk  I  =  v  (kj)«Ak  <  |v  (ki )  I • I Ak | . 

(2.6-6) 

Therefore  three  phonon  processes  are  only  possible  if 
for  some  thermal  phonons  k j , 


v  (k  j )  >  I  v  I  , 

-g  -  -  1 -p 1  sound, 


(2.6-7) 


i.e.  there  must  be  anomalous  dispersion  in  order  that  a 
longitudinal  sound  wave  may  interact  via  three-phonon 
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processes . 

Now  by  intermediate  temperatures  I  mean: 


“V1  and  VTu  • 


(2.6-8) 


This  region  may  not  exist  at  sufficiently  high  frequencies, 
but  it  could  be  extended  if  anomalous  dispersion  is  pre¬ 
sent.  If  (2.6-8)  is  satisfied,  Maris  (1971)  gives: 
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2  o 


i  .  2  2  ' 

1  +  00  T 

0 


where  k  is  the  thermal  conductivity  and 


(2.6-9) 


(2.6-10) 
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(2 . 6-11) 


Several  other  authors  (Guyer  (1966)  ,  Niklasson  (1970) , 
Banerjee  and  Pao  (1974))  have  made  calculations  in  this 
region . 

Both  normal  heat  diffusion  and  second  sound 
propagation  contribute  to  the  attenuation  when  oot  >>1. 
For  a  simple  isotropic  material. 
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(2.6-12) 
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where  t  is  the  relaxation  time  found  from  the  thermal 
conductivity  and  thoroughly  discussed  by  Guyer  and 
Krumhansl  (1966)  and  Hogan  et  al.  (1969).  Figure  4  shows 
the  relaxation  times  x^  ,  x^  and  x  in  helium  at  124.4 
bar  and  for  a  1  cm  wide  specimen.  The  curves  are  extra¬ 
polated  from  the  data  of  Hogan  et  al.  (1969)  and  Lawson  et 
al.  (1973)  in  appendix  2.  Although  for  the  frequencies 
used  in  this  work  we  are  never  in  the  intermediate  temper¬ 
ature  region,  the  relaxation  times  are  not  definitive 
since  different  methods  of  measurement  consistently  give 
different  results  (Rogers  (1972),  Berman  et  al.  (1966), 
Bertman  et  al.  (1966) ) . 

2 . 7  Sound  Propagation  in  the  Low  Temperature  Limit 

In  the  low  temperature  limit  (cox>>1)  we  might 
expect  isothermal  sound  propagation.  However  in  order  to 
even  speak  of  a  sound  wave  we  must  have  a  wavelength  much 
greater  than  the  mean  free  path: 

X  >>  i  or  cox  <<  1.  (2.7-1) 

The  isothermal  region  is  never  reached  since  normal  sound 
propagation  breaks  down  at  the  same  point  that  second 
sound  propagation  breaks  down  (equation  (2.4-24)). 

As  the  relaxation  time  becomes  very  large  the 
phonons  no  longer  interact  with  the  lattice  except  at 


. 


T  ( K ) 


Figure  4 
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boundaries  and  so  they  propagate  ballistically .  The 
sound  wave  may  then  be  treated  as  a  beam  of  low-energy 
phonons  which  populate  a  particular  state  N(KJ).  These 
phonons  are  scattered  and  the  energy  levels  change  only 
through  the  anharmonicity  of  the  crystal.  Maris  (1967, 
1971)  has  done  a  perturbation  calculation  and  finds  that 
this  so-called  zero  sound  velocity  is  given  by 


v 

z 


V  + 

a 


Y  2 
Yeff 


C..T 

v 


2  pVv 


(2.7-2) 


where  .  _ 

2  /Y  (kj)K-v  (kj)  \  2  2 

yeff  =\  tn  -  K-v  (k.nD)+  <Y  <t3)>  -  <T(kj)>  , 

\  -  -g  -d  P/ 

(2.7-3) 

v  is  the  adiabatic  velocity  from  (2.5-6)  and  P 
signifies  principal  part. 

Therefore  there  should  be  a  shift  in  velocity 
near  wt  ^  1.  Together  the  second  and  third  terms  of 
(2.7-3)  are  positive  while  the  first  term  is  positive  for 
normal  dispersion  and  either  positive  or  negative  for 
anomalous  dispersion. 

In  order  to  calculate  attenuation  in  the  low 
temperature  regime,  the  effects  of  scattering  from 
phonons  of  various  polarizations  must  be  considered.  The 
calculations  are  very  lengthy  but  in  general  (Maris  (1971) )  , 


4 

a  ^  S3  T  . 


(2.7-4) 
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2 . 8  Dislocations 

Because  of  solid  helium's  high  purity,  disloca¬ 
tions  and  vacancies  (section  2.9)  are  the  only  lattice 
defects  we  need  consider. 

If  dislocations  are  pinned  at  nodes  (intersections) 
then  they  can  only  slip  between  pinning  points  so  that 
under  the  influence  of  a  low  amplitude  strain  field  we 
have  a  situation  not  unlike  a  vibrating  string.  This 
model  has  been  investigated  by  Granato  and  Lucke  (1956, 
1966)  who  find, 


a  (oo)  = 


K 


2  , 
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2ttv*  r  ,  2  2,2  .  .2.  ' 
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Here  co  =  n /£  , 
o 


n2  =  2G/p(l-v)  =  2vq2/(1-v), 


and  K  =  AA  n  R2 

0 


(2.8-1) 

(2.8-2) 

(2.8-3) 

(2.8-4) 

(2.8-5) 


where  Z  is 
ing  points, 
is  Poisson's 


the  length  of  the  dislocation  between  pinn- 
G  -  c^  is  the  shear  modulus,  v  -  cg3/c33 


,n 


ratio,  d  a  T  ,  where  n  ^  5  (Brailsford 


(1972))  is  a  damping  factor,  A  is  the  total  length  of 
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2  2 

dislocations  per  unit  volume,  Aq  -  4(1-v)/tt  and  R 
is  an  orientation  factor  given  by 


R  =  siny  cosy  cost))  . 


(2.8-6) 


<t>  is  the  angle  between  the  slip  plane  and  slip  direction. 

Not  all  dislocation  loops  are  of  the  same  length  £ 
so  some  distribution  must  be  used  such  as  that  proposed  by 
Koehler  (1952) : 
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(2.8-7) 


(2.8-8) 


(2.8-9) 


Wanner  et  al. (1976)  have  approximated  equation 
(2.8-9)  by  realizing  that  there  are  only  small  contribu¬ 
tions  to  Av/v  from  very  long  loops  (£>>L)  and  very  short 

loops  (£<<£  )  where, 
c 

£  =  n/w  =  [2G/p (l-v)m2]^  .  (2.8-10) 

c 

As  well,  loops  just  slightly  larger  and  slightly  smaller 

than  £  effectively  cancel  one  another.  Therefore  only 
c 
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two  short  ranges  of  £,  one  larger  and  one  smaller  than 

£  ,  have  much  effect  and  we  may  write: 
c 


2  2 
4v  RZ 
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2  .2 
m  +d 
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2  2  2  2 
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(2.8-11) 


In  (2.8-10)  >  A^  and  k  is  an  additional  parameter 

fixed  by  the  dislocation  length  distribution: 


k  = 


£  /<£.  > 

c  1 
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(2 . 8-12) 


For  an  exponential  distribution, 
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(2.8-13) 


(2.8-14) 


and  A-L  =  A  -  A£ 


(2.8-15) 


2 . 9  Vacancies  and  Superfluidity  in  Solid  Helium 

A  recent  review  of  vacancies  in  solid  helium  may 
be  found  in  Glyde  (1976) .  As  originally  shown  by  Hether- 
ington  (1968)  vacancies  in  solid  helium  are  non-localized 
as  a  consequence  of  the  very  large  atomic  motions.  There 
is  then  a  band  of  vacancy  energies  with  its  mean  value 


- 
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being  the  usual  local  vacancy  formation  energy  e  .  The 
total  tunnelling  rate  to  all  other  sites  may  then  be 
estimated  through  the  uncertainty  principle  to  be, 


where  Ae  is  one  half  of  the  bandwidth.  Hetherington 

11  -1  3 

calculated  Ae  ^  6K  and  t  ^  3x10  sec  in  He  .  Values 

4 

for  He  are  similar  (Miyoshi  et  al.  (1970) ) . 

3 

In  He  ,  the  non-localized  vacancies  will  have  a 
very  short  mean  free  path  ^  a,  since  the  uncorrelated 
spins  do  not  provide  a  periodic  background  (Glyde  (1976))  . 

4 

However  in  He  ,  the  rapidly  tunnelling  vacancies  can 
take  on  wave-like  properties  (Andreev  and  Lifshitz  (1969) , 
Guyer  et  al. (1971) ) ,  but  only  when  they  can  travel  several 
lattice  sites  during  one  phonon  period. 


i.e.  T  <  ht/k  =  2Ae/k  . 


(2.9-2) 


Andreev  and  Lifshitz  (1969)  have  estimated  that  vacancies 
will  propagate  as  waves  when  the  diffusion  coefficient  is 
greater  for  non-localized  vacancies  than  for  localized 
ones , 


1/9 


i  .  e . 


T  <  0  ( Ae/0) 


(2.9-3) 
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Also  the  repulsive  interaction  energy  between  vacancies 
must  be  less  than  the  bandwidth  so 

T  <  eQ/in  (mc2/Ae)  ,  (2.9-4) 

where  m  is  the  atomic  mass  and  c  is  the  sound  velocity. 
Then  for  nearest  neighbours  the  energy  spectrum  is 

e  (q)  =  +  2 Trhz t  cos  qa/2  ,  (2.9-5) 

where  zt  again  is  the  bandwidth.  The  calculated  band¬ 
width  is  quite  insensitive  to  volume  (Guyer  et  al.  (1971) ) . 

The  existence  of  the  vacancy  wave  energy  band 
and  of  bound  vacancy  states  (Andreev  (1975))  will  modify 
the  phonon  spectrum  near  the  band  edges  and  at  ±  4Ae 
and  ±  2As  (Kondratenko  (1976)).  However  it  would  be 
fortuitous  if  these  energies  occurred  near  the  sound  wave 
frequency. 

We  might  expect  that  the  large  zero  point  motion 
in  solid  helium  would  result  in  the  presence  of  a  finite 
number  of  vacancies  even  at  T  =  0 ,  so  that  we  may  have  a 
Bose  condensate  and  a  superfluid  solid  (Andreev  and 
Lifshitz  (1969),  Chester  (1970)).  Then  there  will  be  four 
sound  propagation  branches  with  the  three  normal  modes 
having  spectra: 


■ 
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co(k)  =  03q  (k)  +  —  (p^  cos  y'  +  sin  y') 

(2.9-6) 

for  a  hexagonal  crystal.  y'  is  the  angle  between  the 

c-axis  and  the  direction  of  polarization  (y  -  6),  while 
s 

p^  is  the  superfluid  density  in  direction  i: 

P®  =  P  -  [ 92e/9p2] ~1  .  (2.9-7) 

__  g 

Estimates  of  p  /p  range  all  the  way  from  10  (Guyer 

s 

(1971))  to  0.2  (Saslow  (1976)). 

The  existence  of  zero  point  defects  may  not  be 
necessary  (Leggett  (1970) )  for  superfluidity  in  the  solid, 
although  it  would  enhance  the  effect.  Saslow  (1977)  has 
done  an  extensive  microscopic  and  hydrodynamic  treatment 
of  the  problem  from  this  viewpoint.  He  shows  that  a 
superfluid  velocity  cannot  be  developed  from  a  Galilean 
transformation  on  a  periodic  material.  Therefore  it 
would  be  extremely  difficult  to  detect  true  superfluid 
flow.  However  there  will  still  be  a  small  effect  on  the 
sound  velocity: 

v2  =  v02  [1-6  v2/(vo2-v2)]  (2.9-8) 

where  v^  is  the  fourth  sound  velocity  given  by 
v42  =  pS  Oij/3p)s,e 


(2.9-9) 
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and 

y  is 
S  is 


6  =  1  -  p  1  S  ( 9  p/3S )  +  p  1  (3p/3e  ) 

1-i/E  x  x  o  ,  y 

(2.9-10) 

the  chemical  potential,  e  is  the  strain  and 
the  entropy. 


CHAPTER  3 


EXPERIMENT 


3 . 1  Cryogenic  Systems 

The  two  major  experimental  requirements  for  the 
study  of  solid  helium  are  a  refrigeration  system  and  a 
pressure  handling  system. 

3.1-1  Cryostat 

The  low  temperature  end  of  the  cryostat  is  shown 
in  figure  5.  The  pressure  cell  was  bolted  onto  a  graphite 

4 

and  copper  support  column  which  was  attached  to  the  He 
pot.  The  weight  of  the  pressure  cell  did  not  allow  it  to 
be  supported  by  the  delicate  dilution  refrigerator. 

Copper  rings  at  several  stages  provided  convenient  loca¬ 
tions  for  thermal  anchoring. 

Cooling  to  about  1  K  was  obtained  by  circulating 

gas  through  the  dilution  refrigerator  while  pumping  on  the 

4  4 

He  pot.  Under  normal  operating  conditions  the  He 

would  last  for  several  days  without  a  refill.  In  order 

to  reach  lower  temperatures  the  gas  mixture  of  approxi- 

3  4 

mately  0.2  moles  He  and  0.5  moles  He  was  condensed 
through  the  still.  Then  the  still  was  pumped  slowly  until 
the  temperature  of  about  0.5  K  was  reached.  At  this  point 
circulation  was  started,  the  still  heater  turned  on  to 
about  1  mW  and  the  dilution  refrigerator  commenced 
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FIGURE  5 


CRYOSTAT 


1.  Needle  valve 

2.  Pressure  line 

3.  Allen-Bradley  carbon  thermometers  (10  ft, 

4.  Thermal  contacts 

5.  He4  pot 

6.  Still  heater  (200  ft) 

7.  Still  (S.H.E.  minifridge) 

8.  Copper  rings 

9.  Mixing  chamber 

10.  Speer  carbon  thermometers  (220  ft,  0 . 5W) 

11.  Pressure  line  heater  (100  ft) 

12.  Vacuum  pumping  line 

13.  He  pumping  line 

14.  He  pumping  line 

15.  He3  return  line 

16.  Flow-limiting  capillary 

17.  Continuous  heat  exchanger 

18.  Calibrated  germanium  thermometer  (S.H.E 

19.  BeCu  pressure  cell 

20.  Heater  (100  ft)  for  crystal  growth 

21.  rf  feedthrough 


0.1W) 


Corp . ) 


22.  Vacuum  can 
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operation.  When  the  temperature  reached  a  stable  minimum 
the  still  heater  power  was  reduced  to  0.25  mW  and  recon- 
densation  of  He  was  stopped.  The  refrigerator  then 
continued  in  a  one-shot  mode  to  the  lowest  temperatures 
attainable. 

Without  the  presence  of  the  pressure  cell,  its 
support  column  and  the  pressure  line  the  refrigerator 
reached  a  minimum  temperature  of  approximately  40  mK. 
With  the  experiment  in  place,  240  mK  was  the  lowest  pos¬ 
sible  temperature. 

The  cooling  power  of  a  dilution  refrigerator  is 
given  by  Betts  (1968)  as 

0  =  10n3RTD2  ,  (3.1-1) 

where  n^  is  the  circulation  rate  in  moles/sec,  R  is  the 
gas  constant  and  TQ  is  the  temperature  of  the  mixing 
chamber.  The  circulation  rate  was  measured  with  a  flow¬ 
meter  of  the  type  described  by  Almond  et  al  (1972)  to  be 
-5 

6x10  moles/sec.  Then  equation  (3.1-1)  becomes 

Q  =  5xl04  Td2  erg/sec.  (3.1-2) 

Most  of  the  heat  input  was  determined  to  be  flowing 
through  the  solidified  helium  in  the  pressure  line.  From 
Hogan  et  al.  (19  69)  the  thermal  conductivity  for  T  <;  1  K 
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and  in  a  tube  of  inner  diameter  d  is  given  by 

K  =  I  Cv  VD  a  d/V  '  (3.1-3) 

where  a  is  a  numerical  constant  slightly  larger  than  1. 
Then  for  a  line  of  length  L  anchored  at  the  He4  pot  of 
temperature  Tq,  the  heat  flow  is 

Q  -  4.3x10  d3  (Tq4-Td4)/L  erg/sec.  (3.1-4) 

Equating  (3.1-2)  and  (3.1-4),  we  find: 

Td2  =  8 . 6xl02  d3  To4/L  .  (3.1-5) 

For  L  =  40  cm,  d  =  0.03  cm  and  Tq  =  1.2  K,  then  TD  *  80  mK. 
With  all  other  heat  losses  included,  the  experiment  should 
still  cool  to  at  least  100  mK.  The  cause  of  this  large 
discrepancy  was  only  found  after  the  present  data  had  been 
recorded.  Through  a  miscalculation  the  thermal  anchor 

4 

at  the  He  had  far  too  small  a  conductivity  (approximately 
40* AT  erg/sec).  This  is  negligible  compared  to  the  heat 
input  of  2000  erg/sec  calculated  from  (3.1-2)  with 
Td  ^  0.2  K. 

The  illustration  of  the  pumping  systems  in  figure 
6  is  self-explanatory.  I  should  only  add  that  the  rather 
large  volume  (2  litres)  of  the  first  liquid  nitrogen  trap 
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FIGURE  6 


PUMPING  SYSTEMS 

1.  He3-He4  gas  mixture  storage 

2.  Molecular  sieve  (Linde  13X) 

3.  Flowmeter 

4.  Dilution  refrigerator  (S.H.E.  minifridge) 

5.  Booster  pump  (NRC  123) 

3 

6.  Sealed  rotary  pump  for  He  circulation 

3 

7.  He  supply 

4 

8.  He  supply 

9.  Diffusion  pump 

10.  Rotary  pump  for  cleaning  system 

4 

11.  He  exchange  gas  storage 

12.  Diffusion  pump 

13.  Rotary  pump  for  vacuum  can 

14.  Booster  pump  (NRC  123) 

4 

15.  Rotary  pump  for  He  pot 

16.  Thermocouple  gauges  (NRC  531) 

17.  Cold  cathode  gauges  (Balzers  HV5 ) 
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accounted  for  the  quantity  of  gas  required  in  the  circu¬ 
lation  system.  The  volume  of  the  trap  containing  the 
molecular  sieve  was  250  m£,  and  of  the  storage  volume, 

90  £. 

3.1-2  Pressure  Cell 

The  pressure  cell  or  bomb  (figure  7)  was  construc¬ 
ted  of  hardened  beryllium  copper  to  hold  pressures  of 
300  bar  and  to  maintain  the  helium  at  constant  volume. 

The  inner  diameter  is  28  mm  while  the  distance  from  the 
^-inch  transducer  to  the  mirror  can  be  varied  from  7  to  10 
mm.  With  this  arrangement  even  signals  with  the  maximum 
beam  deviation  (see  section  2.3)  can  be  observed  without 
any  losses  through  hitting  the  side  wall.  Such  a  situa¬ 
tion  is  essential  in  order  to  make  quantitative  measure¬ 
ments  of  ultrasonic  attenuation. 

The  plate  from  which  the  transducer  was  suspended 
was  tightly  held  between  a  ledge  and  a  split  ring  circlip 
by  12  set  screws.  Six  screws,  three  pushing  and  three 
PuHin<?'  were  used  to  hang  the  transducer  assembly  from 
this  plate.  The  reflecting  mirror  was  made  of  BeCu  and 
polished  flat  to  h  wavelength  of  sodium  light.  After  its 
threads  were  coated  with  Stycast  2850  GT  epoxy  (Emerson 
and  Cummings  Inc.)  it  was  screwed  into  place.  The  entire 
cell  was  then  bolted  onto  a  copper  plate  on  the  support 
structure  (figure  5),  with  the  transducer  end  down.  The 
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Figure  7  High  pressure  cell  and  ultrasonic  assembly 
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pressure  line  and  rf  miniature  coaxial  cable  (Lake 
Shore  Cryotronics  Type  A)  were  entered  into  the  bomb 
through  separate  high  pressure  feedthroughs. 

The  pressure  capillary  was  a  CuNi  tube  of  0.3  mm 
i.d.  and  0.5  mm  o.d.  leading  as  far  as  the  helium  bath 
and  from  there  a  large  tube  was  led  out  of  the  cryostat  to 
an  external  ballast  volume  (about  2x10^  times  the  volume 
of  the  bomb) .  The  helium  used  in  the  pressure  system  was 
purified  by  passing  it  through  liquid  helium. 

3 . 1-3  Crystal  Growth 

Using  the  following  procedure,  high  quality  crys¬ 
tals  were  obtained  for  about  half  of  the  attempts. 

The  bomb  was  first  warmed  up  to  several  degrees 
above  the  melting  point  and  then  slowly  cooled  by  circu¬ 
lating  gas  through  the  dilution  refrigerator.  The  pressure 
line  was  heated  by  passing  3  ma  through  the  80  ft  manganin 
wire  wound  around  it.  This  heater  had  several  extra  turns 
near  the . bomb  end. 

The  progress  of  the  solidification  process  in  the 
bomb  was  followed  by  monitoring  the  ultrasonic  signal.  As 
solidification  commenced,  the  echoes  visible  in  the  liquid 
disappeared  entirely.  As  the  liquid-solid  interface 
reached  the  transducer,  the  echo  train  for  the  solid 
appeared.  At  this  point  the  gas  circulation  rate  was 
raised  to  increase  the  cooling  rate.  A  few  crystals  were 
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annealed  in  a  temperature  gradient  for  twenty  or  thirty 
minutes  at  this  time  by  turning  on  the  temperature- 
controlled  heater  on  the  warm  end  of  the  bomb.  Any  crys¬ 
tals  with  fewer  than  100  visible  echoes  after  this  pro¬ 
cedure  were  discarded.  With  the  bomb  heater  then  switched 
off,  the  current  through  the  pressure  line  heater  was 
increased  to  about  10  ma  and  solidification  continued. 

This  large  current  was  necessary  to  keep  the  pressure  line 
open  at  all  times  in  order  that  solidification  might  occur 
at  constant  pressure.  Superior  crystals  have  previously 
been  grown  by  this  method  (Franck  and  Hewko  (1973)).  in 
a  few  cases  where  crystals  were  accidentaly  grown  at 
constant  volume,  the  ultrasonic  echoes  almost  entirely 
disappeared.  Once  the  temperature  of  the  bomb  began  to 
fall  rapidly,  solidification  was  assumed  to  be  completed. 
All  heaters  were  then  switched  off  and  gas  circulation  was 
stopped.  The  entire  process  lasted  6-8  hours. 

We  have  no  direct  evidence  for  the  quality  of  the 
crystals,  but  there  are  two  reasons  for  supposing  that  we 
had  generally  high  quality  specimens.  The  large  number 
of  clean  echoes  visible  (as  many  as  175  for  crystal  A8-d) 
imply  a  small  number  of  crystal  defects  which,  if  present, 
would  adversely  affect  the  attenuation.  Some  very  low 
angle  grain  boundaries  could  be  present  but  these  are  of 
little  importance  in  ultrasonics  as  long  as  they  are  not 
too  numerous.  As  well,  the  range  of  absolute  velocities 
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observed  (including  discarded  crystals)  was  735-840  m/sec 

this  spans  almost  the  entire  range  of  velocities  expected 
(see  figure  2 ) . 

3 . 2  Thermometry 

The  principal  secondary  thermometers  were  con¬ 
structed  from  220  ft,  0.5  Speer  carbon  resistors  (Black  et 
al.  (1964) ,  Oda  et  al.  (1974) )  ,  whose  resistance  was 
measured  by  an  S.H.E.  model  ARB  automatic  resistance 
bridge.  Power  dissipation  was  typically  10-11  -  10~12 
watts.  The  thermometers  were  prepared  by  grinding  the 
outer  coatings  from  the  resistors  and  covering  them  with 
thinned  GE  varnish  (#7031)  which  was  baked  on  at  90°C. 

The  resistors  were  inserted  and  varnished  into  machined 
copper  blocks  which  were  screwed  to  the  copper  support 
plate,  to  the  bomb  and  to  the  still.  The  42  AWG  manganin 
wire  leads  were  thermally  anchored  at  the  mixing  chamber, 
the  heat  exchangers,  the  still  and  the  He4  pot. 

The  thermometers  were  calibrated  in  two  ways. 
Initially  a  CMN  paramagnetic  thermometer  was  constructed 
(Abel  et  al.  (1964)  ,  Fisher  (1972)  )  and  the  susceptibility 
measured  with  a  Hartshorn-type  mutual  inductance  bridge 
(Maxwell  (1965)).  The  CMN  salt  pill  was  calibrated 
against  the  vapour  pressure  of  He4  at  temperatures  above 
IK.  As  there  were  some  problems  reading  the  bridge  this 
method  was  not  entirely  satisfactory  so  the  thermometers 
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were  later  compared  to  an  S.H.E.  germanium  resistance 
thermometer  calibrated  from  30  mK  to  5  K  and  measured  on 
the  ARB  bridge.  This  thermometer  was  held  in  a  machined 
copper  block  with  Corning  silicone  vacuum  grease.  This 
later  calibration  raises  the  temperatures  recorded  with 
the  data  in  this  work  by  about  40  mK. 

3.3  Ultrasonics 


3 . 3-1  Transducer  Assembly 

The  ultrasonic  etalon,  as  supplied  by  Valpey- 
Fisher,  was  a  ^-inch  diameter  quartz  transducer  with  a 
fundamental  frequency  of  5  MHz.  It  was  overtone  polished 
to  provide  good  excitation  of  higher  harmonics,  and  coated 
on  both  sides  with  a  thin  film  of  AuCr  alloy.  As  can  be 
seen  in  figure  7,  the  front  of  the  transducer  was  held  at 
the  rim  by  a  brass  flange  3  mm  thick  and  25  mm  in  diameter 
with  a  10  mm  diameter  hole  in  the  centre.  Six  radial 
grooves  were  scored  on  the  back  of  the  flange  to  prevent 
an  unwanted  pressure  build-up.  The  flange  and  the  bottom 
end  of  the  main  body  of  the  transducer  assembly  were  pol¬ 
ished  flat  to  2  fringes  of  sodium  light.  Six  screws  held 
the  flange  to  the  body.  Backing  the  transducer  was  a 
brass  piston  polished  to  2  fringes  of  sodium  light  and 
inserted  into  a  teflon  sleeve.  The  piston  was  bevelled 
at  the  top  to  dampen  reflections  in  the  brass  while  the 
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teflon  sleeve  had  a  2°  bevel  so  that  it  would  be  self¬ 
levelling.  With  the  teflon  sleeve  in  place  the  piston 
had  the  rf  lead  wire  attached  to  it  and  then  was  spring- 
loaded  into  the  large  brass  cylinder  (25  iron  diameter) 
which  formed  the  main  body  of  the  assembly. 

Optically  polished  surfaces  and  the  choice  of  a 
spring  of  the  proper  size  and  force  constant  were  found 
to  be  crucial  for  the  attainment  of  a  good  ultrasonic 
signal. 

The  assembly  was  mounted  within  the  bomb  by  three 
pairs  of  push-pull  screws  as  described  in  section  3.1-2. 

3 . 3-2  Electronics 

Figure  8  is  a  block  diagram  of  the  electronics. 
The  technique  employed  was  a  modified  pulse-echo-overlap 
method  as  described  by  May  (1958),  Chung  et  al.  (1969), 
Papadakis  (1964,  1967)  and  Van  Nest  et  al.  (1969).  The 
heart  of  the  system  is  a  highly  stable  (to  better  than 

7 

1  part  in  10  per  day)  frequency  synthesizer,  model  5100 
by  Rockland  Instruments.  It  operates  at  a  frequency  v 
in  the  neighbourhood  of  50  kHz  such  that 

1/v  =  2D/v  ,  (3.3-1) 

where  v  is  the  sound  velocity  and  D  is  the  distance 
from  transducer  to  mirror.  This  signal  is  divided  by  a 
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BLOCK  DIAGRAM  OF  ULTRASONICS 


Figure  8 
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factor  of  100  or  1000  in  a  MATEC  121  decade  divider  - 
dual  delay  generator  and  is  used  to  trigger  the  rf  pulse 
generator-wide  band  amplifier  (MATEC  6600  with  plug-in 
760)  .  The  rf  pulse  is  sent  by  50  Q,  coaxial  cable  through 
a  system  of  two  passive  tuning  networks,  one  supplied  by 
MATEC  (model  60)  and  the  second  a  twin-tee  stub  tuner 
network  (Wanner  (1970) ) .  The  pulse  is  sent  to  the  trans¬ 
ducer  and  the  returning  echoes  again  pass  through  the 
tuning  networks  to  be  amplified  by  the  model  6600.  The 
resulting  signal  is  applied  to  the  y-axis  of  a  Hewlett- 
Packard  182C  oscilloscope  with  plug-ins  1808A  and  1825A. 
The  oscilloscope  is  triggered  with  the  undivided  pulse¬ 
shaped  signal  from  the  frequency  synthesizer.  If  one  uses 
the  continuous  delay  on  the  oscilloscope  trigger  circuit, 
any  part  of  the  pulse  train  may  be  observed.  Because  of 
the  image  retention  of  the  oscilloscope  screen  and  the 
human  eye,  the  rf  contents  of  all  echoes  may  be  observed 
simultaneously  and  will  all  appear  to  be  lined  up  if 
equation  (3.3-1)  is  satisfied  (as  in  figure  9(b)  ).  Then 

the  velocity  is  2Dv.  If  the  velocity  or  trigger  frequency 
changes  then  the  echoes  shift  relative  to  one  another  as 
seen  in  figure  9(c,d).  As  well,  the  oscilloscope  may  be 
internally  triggered  in  order  to  view  the  entire  echo 
train  (figure  9(a)). 

When  the  signal-to-noise  ratio  became  small,  a 
double  negative  strobe  was  applied  to  the  z-axis  of  the 


. 
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Figure  9  -  Ultrasonic  Echo  Train  and  Echo  Overlap  (in  Liquid  Helium) 
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oscilloscope.  The  two  strobe  pulses  were  delayed  by  the 
dual  delay  generator  so  that  they  were  synchronous  with 
any  two  selected  echoes.  Then  only  those  two  echoes 
remained  visible  on  the  oscilloscope  screen  while  other 
echoes  and  most  of  the  noise  were  deleted. 

Some  attempt  was  made  to  measure  attenuation  with 
a  MATEC  2470A  attenuation  recorder.  It  automatically 
measures  the  decrement  between  two  pre-selected  echoes  and 
continuously  copies  it  onto  a  chart  recorder.  The  presence 
of  the  highly  non-exponential  echo  envelope  (see  section 
3.3-4)  limited  the  usefulness  of  this  method.  Instead  the 
attenuation  was  determined  by  photographing  the  echo  train 
with  a  Hewlett-Packard  197A  oscilloscope  camera  and  manu¬ 
ally  measuring  the  decrement. 

The  initial  pulse  had  an  amplitude  of  200  to  400 
Vpp  and  a  width  of  approximately  4  ys.  The  repetition 
rate  was  v/100  at  higher  temperatures  and  v/1000  at  lower 
temperatures.  There  was  no  discernible  power  dependence 
of  the  velocity  upon  changing  either  the  pulse  height  or 
repetition  rate. 

Very  careful  tuning  of  both  passive  networks  was 
necessary  in  order  to  obtain  the  best  signal  and  partially 
compensate  for  the  very  large  acoustic  mismatch  (over  100) 
between  quartz  and  solid  helium. 
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3.3- 3  Method 

Initially  the  cell  was  filled  with  methanol  and 
the  resulting  ultrasonic  echoes  were  observed.  The  de¬ 
gree  of  non-parallelism  between  the  transducer  and  reflec¬ 
tor  was  monitored  by  watching  the  echo  envelope  (section 

3.3- 4) .  The  push-pull  screws  were  then  adjusted  and 
tightened  to  lock  the  transducer  in  place,  parallel  to  the 
mirror . 

Measurements  were  taken  first  in  liquid  helium 

while  the  cell  was  maintained  in  equilibrium  with  the 

helium  bath  through  the  medium  of  He^  exchange  gas  in  the 

vacuum  can.  From  the  known  velocity  v  in  the  liquid 

(Vignos  and  Fairbank  (1966)  as  corrected  by  Wanner  (1970) 

and  Abraham  et  al. (1970) ) ,  the  transducer  to  reflector 

distance  D  was  determined  from  D  =  v  /2v.  It  was  assumed 

Li 

that  D  changed  little  upon  freezing  so  that  the  absolute 
velocity  could  then  be  determined  in  the  solid.  Thence¬ 
forth  only  changes  in  velocity,  or  changes  in  v,  needed 
to  be  recorded. 

The  true  condition  for  overlap  is  given  by 
McSkimin  and  Andreatch  (1962)  as. 


1  _  2pD  n 
v  v  f 


(3.3-2) 


Here  p  is  the  number  of  round  trips  between  trigger 
pulses  (usually  p  =  1) ,  f  is  the  rf  frequency  and  n  =  0 
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when  succeeding  echoes  are  correctly  superposed.  The 
oscilloscope  can  be  triggered  with  n  ^  0  so  that  each 
echo  is  displaced  by  n  rf  periods  relative  to  the  pre¬ 
vious  echo.  The  rf  oscillations  will  still  appear  aligned 
and  for  long  pulses  or  short  echo  trains  the  effect  will 
not  be  obvious.  However  from  (3.3—2)  it  can  be  seen  that 
if  f  is  varied,  the  overlap  condition  will  not  be  ful- 
ft  7^  McSkimin '  s  method  then  is  to  vary  f 

for  several  overlaps  in  order  to  find  the  one  independent 
of  f.  in  practice  this  choice  was  always  confirmed  by 
comparing  1/v  with  the  average  time  between  echoes  along 
the  entire  pulse  train.  McSkimin' s  criterion  was  applied 
at  each  different  rf  frequency. 

Once  the  absolute  velocity  was  determined,  the 
temperature  was  slowly  lowered  and  v  was  adjusted  regu¬ 
larly  in  order  to  keep  the  echoes  in  alignment  and  to 
follow  the  same  rf  oscillations  over  the  entire  temperature 
range.  Measurements  for  any  one  rf  frequency  were  per¬ 
formed  over  a  period  of  two  days,  often  repetitively  over 
certain  temperature  ranges.  Some  crystals  were  cycled 
with  respect  to  temperature  and  frequency  with  excellent 
reproducibility . 

3 . 3-4  Errors  in  Measurement 

The  uncertainty  in  the  measurement  of  D,  and 
therefore  in  the  absolute  velocity,  is  0.1%,  based  on  the 
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reproducibility  of  the  measurements.  The  largest  errors 
in  the  determination  of  the  orientation  arise  from  the 
extrapolation  of  the  elastic  constants  to  17.4  cm3 
(appendix  1).  This  is  estimated  to  be  about  2%,  since 
the  extrapolation  to  16  cm3  (Reese  et  al.  (1971))  is  correct 
to  better  than  5  %.  The  error  in  orientation  is  very 
much  a  function  of  the  velocity,  because  of  the  wide  vari¬ 
ation  in  dv/dy  (figure  2).  For  42°  <  y  <  90°,  the  velocity 
is  double-valued.  In  these  cases  a  most  likely  orientation 
is  chosen  by  observing  the  interference  pattern.  A  pattern 
similar  to  that  seen  in  the  liquid  was  considered  to  arise 
from  an  orientation  with  a  small  beam  deviation  (figure  3) . 
For  orientations  with  a  large  A,  the  beam  will  sample  a 
different  part  of  the  reflector  and  so  the  interference 
pattern  will  change.  The  measuring  error  in  Av/v  varies 
inversely  with  the  number  of  echoes  observed.  In  practice 
the  error  for  100  echoes  is  about  2xl0~6  and  for  2  echoes 

_  c 

(as  measured  at  low  temperatures)  about  5x10  .  The 

limiting  factors  are  the  jitter  in  the  trigger  pulses  and 
various  small  effects  which  are  dependent  upon  echo 
number. 

The  finite  pulse  length  causes  variations  in 

transit  time  for  different  parts  of  the  pulse  (Merkulov 

and  Tret'yakov  (1975))  so  that  there  was  an  observed 
•  -5 

velocity  change  of  about  2-3x10  across  the  pulse.  This 
error  is  effectively  eliminated  by  always  following  the 


63 


same  rf  oscillations  near  the  centre  of  the  pulse  at  all 
temperatures . 

Also  exerting  a  small  influence  on  the  sound 
propagation  is  diffraction  of  the  beam.  This  effect  has 
been  investigated  by  a  number  of  authors  (Seki  et  al 
(1956),  Khimunin  (1975),  Papadakis  (1966,  1972)).  Papada- 
kis  (1966)  has  made  calculations  for  anisotropic  materials 
so  his  results  are  of  the  greatest  interest  to  the  present 
work.  Diffraction  of  the  beam  results  in  a  wave  front 
which  is  no  longer  planar,  and  varying  with  distance 
travelled,  so  that  there  is  a  phase  shift  in  the  detected 
signal.  As  well  the  wave  will  interfere  with  itself  at 
the  transducer,  causing  an  apparent  change  in  attenuation. 
This  last  effect  is  negligible  when  compared  with  other 
problems  involved  in  measuring  the  attenuation  in  solid 
helium.  In  general  the  phase  shift  is  a  complicated 
function  of  wavelength  A,  transducer  diameter  a  and 
distance  travelled  z,  but  the  phase  shift  is  approxi¬ 
mately  linear  in  S  =  zA/a2  for  S  <  0.5  (at  5  MHz, 

S  -  0.023  n,  where  n  is  the  number  of  round  trips) .  The 
proportionality  constant  varies  with  the  anisotropy,  but 
is  ailways  smaller  in  solid  helium  (where  the  anisotropy 

parameter  b  =  v3  of  Waterman  (1959)  is  0.25  along  the 

2 

c-axis  and  0.12  cos  \p  in  the  basal  plane,  ^  being  the 
azimuthal  angle)  than  in  isotropic  liquid  helium  (b  =  0) . 
The  largest  influence  of  diffraction  in  the  present  work 
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is  on  the  measurements  of  dispersion.  The  changes  in 
Av/v  in  the  liquid  should  be  about  -4x10  ^  between  5  and  15 
and  about  -6x10  between  15  and  25  MHz.  The  effect  in 
the  solid  will  always  be  less  and  rarely  more  than  half 
of  that  in  the  liquid. 

The  greatest  geometrical  effect,  especially  with 
respect  to  attenuation  measurements,  arises  from  a  lack 
of  parallelism  between  the  transducer  and  the  reflector. 
This  problem  has  been  investigated  in  a  rather  simple  way 
by  Truell  et  al.  (1969)  and  Abraham  et  al.  (1969)  .  The 
present  calculations  are  much  more  detailed. 

Rigorously  we  begin  with  an  initial  velocity 
distribution  on  the  transducer  Im  f  (x,y;t) .  Then  the 
Fourier  transform  is: 


Fq  ( r ,  0  ;  03 )  = 


e  iwtf  (r,0;t)dt. 
o 


(3.3-3) 


The  phase*  delay  of  the  n —  echo  for  attenuation  a  and 
velocity  c  is: 


^  (r,  0  ;  03)  -  (w/c  “  iot)  x  (distance  travelled)  . 

(3.3-4) 

Then  the  Fourier  transform  for  the  n—  echo  is: 


Fn(r,0;co)  =  FQ(r,0;oo)  e  i(^n  '  0  ;  ^  . 


(3.3-5) 
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Finally,  the 


Im  f  (t)  = 
n 


where  a  is 


received  signal  for  the  n—  echo  will  be: 


[  loot  4 
o—  dwe  — — 

2tt  J  2 

Tra 


a/ 2  2tt 


rdr 


dGF^  (r ,  6  ;  oo )  , 


(3.3-6) 


the  diameter  of  the  transducer. 


Generally,  for  reasonably  long  pulses,  the  centre 
of  the  pulse  may  be  treated  as  monochromatic  so  that  the 
preceding  method  is  only  required  when  considering  the 

ends  of  the  pulse.  Then  we  may  write,  for  a  fundamental 
frequency  ft  , 


a/ 2 


2  TT 


f  (t) 
n  ' 


Tra 


ro.r 


o 


J 

0 


d0e  llf'n(r'9;S,fo(r,e;t)  . 


(3.3-7) 


In  the  simplest  case. 


f  (r , 6  ;  t ) 


0<t<T 

otherwise . 


(3.3-8) 


If  the  reflector  is  tilted  relative  to  the  transducer,  as 
in  figure  10(a),  then: 


<f>  (r,0;o))  =  m(t  +2t*  —  cos0) 
n  n  n  a 


where  t  =  2nD/c 

n  7 


v 

n 


i  ( v  +  2v  '  —  cos 0  )  , 
n  n  a 

(3.3-9) 

2naD 


v ' 

n 


=  nd/c 


nad. 


Transducer  -  Reflector  Arrangements 
(a)  Tilted  Reflector 


< - a - ► 


(b)  Curved  Reflector 


< - a - ► 


Figure  10  Transducer-reflector  configurations 
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Then  equation  (3.3-7)  may  be  integrated  to  give 


-v 


Im  f  (t) 
n 


=  e 


n  2 | J_  (ft  t' ) 
1  1  on 

I  ft  t 1 

1  o  n 


sin[ft(t-t  )  +6-  +  -  ^ 


n 


V  ]* 


(3.3-10) 


Here  6,  =  arctan  co  /ft  -  w  /ft 

1  o  o 


a)  =  ac 
o 


h  (1) 
n 


arg  J  (ft  t 1 ) 
Ion 


ft  = 


ft-iw  , 
o 


and  Jn  is  the  Bessel  function  of  order  n.  Then  the  rf 

-v 

wave  is  modulated  by  an  amplitude  2e  n|j,  (ft  t 1 )  I / I  ft  t'l 

1  1  o  n  1  1  o  n 1 

and  phase-shifted  by  6^  +  0n  ^  .  T^e  maxama  and  minima 
of  this  envelope  may  be  calculated  exactly,  but  since  the 
effects  are  only  important  when  a  is  small,  we  may 
simply  write  that  to  first  order  in  oo^/ft,  the  maxima 
occur  at 


J,  (ftt' )  ftt' 

1  n  _  n 

J  (ftt ' )  v  +2  '  (3.3-11) 

o  n  n 


and  the  minima  at 


J,  (ftf )  =  0. 
1  n 


(3.3-12) 
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if  vn<<2,  then  we  have 


nmax  fid  ^*14,  8.42,  11.62,  14.8,  ...), 

(3.3-13) 

nmin  ~  fid  ^«83,  7.02,  10.17,  13.32,  ...)  . 


Then  for  any  two  selected  echoes  n^  and  n0,  the  true 
attenuation  a  is  given  in  terms  of  the  "measured" 
attenuation  a'  as. 


a  =  a '  + 


£n[J1(t^  )/J1(tj;  )]  +  £n(ni/n2) 
 2  1 


2D (n2-n1) 


(3.3-14) 


To  first  order  in  oj  /fi, 

o 


tan  n 


(1) 


n 


<“0/fl)  &  -  «t^Jo  (Sitp  /Jx  (at')  1 . 


(3.3-15) 


Then  for  two  echoes  n^  and  n2 ,  the  measured  velocity  is 
larger  than  the  actual  velocity  by 


c  2fiD  fi 

The  transducer  does  not  necessarily  vibrate  in  a 
piston-like  manner  as  in  equation  (3.3-8)  (Papadakis 
(1971)),  especially  if  clamped  at  the  edge.  A  simple  case 
to  treat  mathematically  is: 


J  fit' 
o  n. 


J  fit' 

Utn  J  fit'  J~  fit* 

1  1  n..  2  1  n„ 


(3.3-16) 
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fo(r;t)  = 


(l-4r2/a2)  elfit'  0<t<T 
0,  otherwise. 


(3.3-17) 


-v 


n  4  J  (ft  t '  ) 

Then  Im  (t)  =  e  *  - sin  [^(t-t  )+60+n  2  ]  ' 

|fi  t"2  n  2  'n 


n 


o  n 


(3.3-18) 


where  6  =  arctan  [2fito  /  (fi2+u,  2)]  ^200  /Si 

o  0  0 


,  (2) 

and  n  '  arg  J0 (fit')  , 

n  2  0  n 


Then  the  maxima  occur  at 


J2  (fit^/j^fit^) 


fit '/  (v  +4) , 
n  n 


(3.3-19) 


and  the  minima  at 


J0  (fit' )  =  0. 
z  n 


(3.3-20) 


nmin  =  [5.14,  8.42,  11.62,  14.80,  ...] 


(3.3-21) 


and  n 


fi 

max  cd 


[0,  6.38,  9.76,  13.02,  16.22,  ...]. 


Also,  a  =  a'  + 


£n[J2(fie  )/J1(fit^  ) ] +2  £n(ni/n2) 
2  1 

2D  (n2-n1l 


(3.3-22) 
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fitn„  J0 (fit'  ) 
2  2  n  ' 


(3.3-23) 


Equation  (3.3-17)  is  not  necessarily  an  accurate  descrip¬ 
tion  of  the  true  physical  situation,  but  (3.3-18)  suggests 
that  qualitatively  any  deviation  from  a  piston-like  velo¬ 
city  distribution  will  result  in  an  even  larger  apparent 
attenuation  than  that  given  by  (3.3-10). 

The  velocity  corrections  61  and  6^  are  negligible 
when  measuring  dispersion  since  w  <<  ft.  The  error  in 


0 


velocity  between  successive  echoes  (equation  (3.3-16)  or 


(3.3-23))  is  less  than  5xl0~6  for  a  -  0.1  cm  1  and  for 


low  n . 

Instead  of  being  tilted,  the  reflector  may  be 
curved  as  in  figure  10  (b)  where  the  curvature  of  the  reflec¬ 
tor  is  given  by 


(3.3-24) 


Then  the  angle  of  reflection  0  with  the  normal  to  the 

n 


th 

transducer  for  the  n  echo  can  be  shown  to  obey  the 


recursion  relation 


(3.3-25) 


with  0 


o 


0 


71 


and 


0.,  -  26(r^+r^+2rr  cos0)  ^  , 
i  oo 


where  rQ  arises  from  the  beam  deviation  (section  2.3) 


rQ  =  D  tan  A. 


(3.3-26) 


Let  en2  =  0^  un2  (1-2BD)/4BD(1-BD) 


(3.3-27) 


Then  Un  -  2(1-2BD)  -  U  . 


(3.3-28) 


where  U  =  0 


^  =  [1-U-28D)  2]h 


Equation  (3.3-28)  defines  the  Chebyshev  polynomials  of 
the  second  kind.  Then  the  phase  shift  for  the  nfc^ 
reflection  is 


2  2 

,'r  +r  +2rr  cos0 

♦n(r,0;u))  =  -  [2nD  +(  4DU-BD° - 1  x 


n 


x  I  [U  (1-2BD)  -  U  (1-2BD)]2)-  iv  , 
m=l  m  m  1  n 


(3.3-29) 


where  we  ignore  the  slight  variation  in  attenuation  from 
differences  in  reflector  to  transducer  distance.  Then 
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using  (3.3-8)  and  (3.3-29)  in  (3.3-7), 


-v 


n 


f„(t)  =  M  “  ein(t-tn) 


n 


a/2 

dr.r  exp 


o 


{•‘i 


2  2 
r  +r 

o 

4D (l-$D)  x 


J1(VUm-l)  }  J0  _c-  JD(1-°D)J  'Wl1 


(3.3-30) 


If  we  define 


ft 


n 


~  2cD  l  [U  (1-2BD)  -U  l(1-2f?D>J2  ' 
m=l  111  x 


and  use  a  series  representation  for  J  ,  then  it  i 


(3.3-31) 


0 •  -  --  -s  easy  to 


show  (Hansen  (1975)),  that  fn(t)  becomes: 


a  ic  a 

im  f_(t)  =  2|_ "  |P  ( — § — 12,1?  r2) 

4  no 


n 


a2? 


n 


sin  [ft  (t-t  )+n  ^  ] 

n  'n  J 

(3.3-32) 

where  P  is  the  non-central  chi-square  probability 
function  (Abramowitz  et  al.  (1964))  and 


%(3)4~ 

n  2 


arg  P  (-^!_|  2,  i?nr2) 


(3.3-33) 


P  is  given  by 


I 
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P(itna2/4|2,iCnr2)  =  l-e-i?n(ro/2  +  a2/8> 


X 


(i?nr2/2)k(iC„a2/8) j 


I  l 

k=o  j=o 


n 


k!  j  i 


(3.3-34) 


P  is  difficult  to  evaluate,  but  for  r  <<a/2. 


o 


P  -  l-e'i?na  /8  [1  -  ?2r2a2/16], 

no 


(3.3-35) 


Therefore , 


-  2  |  sin  cn  aV16|  [1  - 


2  2 
C  r  £  a 

n  o  ^n 

4  *  8 


(3.3-36) 


Then  the  maxima  of  I P I /r  occur  at: 

n 


tan<;naVl6  =  (?naVl6)  [1  -  ?2r2a2/16  ]  , 


(3.3-37) 


or 


Cn  =  (16/a2)  (l-16r2/a2)  (0,  4.49,  7.73,  10.90,  14.07,  ...) 

(3.3-38) 


Now  r  =  -iL-  l  [U  -U  . ] 
n  2cD  L.  m  m-1 
m=l 


n 


=  ~  [n  +  U2n  (1-3D)/2U12(1-3D)] 
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§  [nB  +  ^  sin  2n/23D  ] 


(3.3-39) 


Therefore  to  find  3,  solve: 


sin  2n/2gD  =  46D  [-n  +  (l-16r2/a2) (0,  4.49,  7.73,  ...)], 


3fia‘ 


where  n  refers  to  a  maximal  echo. 


The  minima  in  |p|/c  occur  when  e  a2/16  =  mir ,  or 

n  n 


(3.3-40) 


sin  2n/2  3D  =  43D  [-n  +  (3.14,  6.28  ,  9.42,  12.57,  ...)]. 


3fia‘ 


Finally,  using  equation  (3.3-36), 


(  3  )  7T 

n  =  77  -  arctan 
n  2 


r  2  1  T  (^roa 

[cot(cna2/16)j  I  1+  -2-§_ 


<?r,roa/4)  2cosCna2/8 


cose  a  /8 
n 


f] 


(3.3-42) 


CHAPTER  4 


RESULTS 


4 . 1  Velocity  of  Sound 

4  ♦  1-1  High  Temperature  Region 

All  crystals  were  measured  at  a  volume;  of  17.4  ± 

3 

0.05  cm  /mole.  The  largest  error  in  the  determination  of 
the  absolute  velocity  arises  from  any  changes  in  the 
transducer  to  reflector  distance  D.  This  error  was 
estimated  from  measurements  of  D  on  successive  crystals 
to  be  0.1%.  The  accuracy  of  the  velocity  for  purposes 
of  orientation  determination  is  estimated  to  be  better 
than  2%. 

Curves  of  velocity  v  (divided  by  an  arbitrary 
v0)  versus  temperature  are  shown  for  ten  crystals  in 
figures  11  to  20.  For  the  first  four  samples,  measure¬ 
ments  were  only  taken  at  5  MHz,  while  all  three  frequen¬ 
cies  of  5,  15  and  25  MHz  were  used  for  the  last  six. 

The  curves  are  least-squares  fits  of  equation 
(2.5-10)  to  the  data  points  in  the  high  temperature  region 
only.  Equation  (2.5-10)  was  obeyed  very  well  in  all  cases 
with  the  constants  having  the  values  listed  in  table  1 
and  illustrated  in  figure  21. 

As  is  apparent  from  figure  21,  the  scatter  in  the 
values  of  a  and  b  is  too  large  to  allow  a  meaningful 


75 


Figure  11  Reduced  sound  velocity  as  a  function  of  temperature. 
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Figure  12  Reduced  sound  velocity  as  a  function  of  temperature. 


Figure  13  Reduced  sound  velocity  as  a  function  of  temperature. 
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Figure  14  Reduced  sound  velocity  as  a  function  of  temperature. 
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Figure  15  Reduced  sound  velocity  as  a  function  of  temperature  and  frequency. 
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Figure  16  Reduced  sound  velocity  as  a  function  of  temperature  and  frequency. 
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Figure  17  Reduced  sound  velocity  as  a  function  of  temperature  and  frequency. 
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Figure  18  Reduced  sound  velocity  as  a  function  of  temperature  and  frequency. 


Figure  19  Reduced  sound  velocity  as  a  function  of  temperature  and  frequency. 
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Figure  20  Reduced  sound  velocity  as  a  function  of  temperature  and  frequency. 
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Orientation  Dependence  of  Velocity 
Coefficients  (5  MHz) 


Figure  21  Orientation  dependence  of  adiabatic  sound  velocity  coefficients  in 
equation  (2.5-10). 
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TABLE  1 

Crystal  y(°)  rf  Freq.  v  (T=0  a  b 

(MHz)  (cm/s)  (cm/s-K  )  (cm/s-IC) 


Al-b 

43 

5 

76416 

-4.35 

0.076 

A2-b 

29 

5 

80951 

-3.45 

0.052 

A4-d 

20 

5 

83447 

-1.74 

-0.123 

A5-a 

24 

5 

82392 

-1.53 

-0.082 

A6-b 

1 

►  46 

5 

75450 

-4.39 

0.076 

A6-c 

J 

I 

15 

75396 

-3.76 

0.013 

A6-e 

> 

>  69 

5 

75112 

-2.91 

-0 . 071 

A6-  f 

J 

I 

15 

75048 

-4 . 42 

0.062 

A7-a 

5 

76176 

-3.48 

-0.067 

A7-b 

►  82 

15 

76160 

-3.93 

-0.017 

A7-c 

25 

76160 

-4 . 54 

0.046 

A7-d 

> 

5 

76027 

-4 . 41 

0.035 

A7-e 

►  42 

15 

75998 

-4.76 

0.072 

A7-f 

J 

25 

75990 

-4 . 02 

-0.003 

A8-a 

5 

75856 

-4 . 00 

-0.003 

A8-b 

>  78 

15 

75827 

-4 . 13 

0.032 

A8-c 

J 

25 

75839 

-4.19 

0.028 

A8-d 

5 

75970 

-3.79 

-0.006 

A8-e 

81 

15 

75948 

-4.27 

0.048 

4-1 

1 

CO 

< 

> 

25 

75963 

-4 . 57 

0.071 

calculation  of  the  temperature  dependence  of  the  elastic 
constants.  However  both  a  and  b  clearly  reach 
maximum  values  at  intermediate  angles  with  minima  at 

y  =  0°  and  90°.  Very  approximately  it  can  be  estimated 
that. 


Ac33/c33  -  -  (3±1)x10~5T4  -  (3±1)x10~6T6  (4.1-1) 

and  &cll/c11  *  -  ( 8±2 ) xlO~5T4  -  (2±l)xlO-6T6  .  (4.1-2) 
^C44/C44  will  have  a  very  large  term,  as  already  found 

3 

in  bcc  He  (Wanner  et  al. (1973) ) .  Table  1  also  shows  a 
fairly  consistent  trend  to  larger  T  contributions  and 

towards  more  positive-going  terms  at  higher  frequencies 

-4  -4 

A  very  large  dispersion  of  -2x10  to  -8x10 
(as  extrapolated  to  T  =  0)  between  5  and  15  MHz,  and 
somewhat  less  between  15  and  25  MHz,  is  observed  over  the 
entire  measured  temperature  range  below  melting.  There 
is  a  slight  tendency  for  the  dispersion  to  become  smaller 
at  higher  temperatures.  In  most  crystals  the  dispersion 
is  normal,  but  because  it  is  some  five  orders  of  magni¬ 
tude  larger  than  expected  from  a  simple  sine  curve,  the 
dispersion  must  necessarily  become  anomalous  at  some 
higher  frequency.  In  fact  two  crystals,  A8-a  and  A8-d 
(figures  19  and  20) ,  do  show  positive  dispersion  at  25  MHz 
To  ensure  that  the  observed  frequency  dependence  was  not 


. 
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an  instrumental  effect,  measurements  were  also  taken  in 

the  liquid  at  4.1  K,  where  the  dispersion  was  found  to 

—  5  _  c 

be  -6x10  between  5  and  15  MHz  and  -1x10  between 

15  and  25  MHz.  These  results  agree  well  with  the  expected 

effects  of  beam  diffraction  (section  3.3). 

The  dispersion  has  no  discernible  dependence 

upon  orientation. 

4.1-2  Low  Temperature  Region 

Although  the  dispersion  continues  to  the  lowest 
temperatures,  in  most  crystals  it  changes  dramatically  as 
the  temperature  is  reduced  because  of  an  anomaly  in  the 
sound  velocity  which  begins  at  some  temperature  between 
1  and  1.8  K  (at  5  MHz)  and  continues  to  the  lowest 
measured  temperature.  The  anomaly  has  been  observed 
previously  (Franck  and  Hewko  (1973),  Wanner  and  Mueller 
(1974),  Wanner  et  al. (1976)).  Within  the  anomalous  region 
the  deviation  from  the  extrapolated  adiabatic  velocity  is 

-3 

as  much  as  Av/v  =  3x10  .  The  anomaly  can  be  either 

positive-going  (the  most  common  situation),  or  negative, 
as  seen  in  the  two  samples  Al-b  and  A8-a  (figures  11  and 
19).  Perhaps  most  striking  of  all  are  crystals  A6-b 
and  A7-a  (figures  15  and  17)  which  show  no  signs  of  any 
anomaly  at  all. 

Neither  the  magnitudes  nor  the  onset  temperatures 
of  the  anomalies  exhibit  any  dependence  upon  orientation. 


. 


Many  crystals  show  a  small  negative-going  anomaly  in  the 
temperature  range  just  above  the  onset  temperature.  In 
most  cases  the  velocity  reaches  a  constant  value  after 
experiencing  a  very  rapid  change.  Crystals  A8-a  and  A8-d 
show  a  velocity  beginning  to  return  towards  the  adiabatic 
value  at  the  lowest  temperature.  These  are  the  same  two 
samples  in  which  anomalous  dispersion  was  observed  at 
25MHz. 

The  anomalies  exhibit  a  marked  frequency  depend¬ 
ence.  At  15  MHz  they  become  smaller  and  move  to  tempera¬ 
tures  some  15  to  35%  lower,  and  at  25  MHz  are  smaller 
again  and  at  the  same  or  a  slightly  lower  temperature. 

The  small  high  temperature  anomaly  becomes  larger  at 
higher  frequencies,  especially  in  crystal  A8-d. 

4.2  Attenuation  of  Sound 


4.2-1  The  Echo  Envelope 

In  order  to  obtain  a  reliable  value  for  the 
attenuation  when  it  is  small,  the  effects  of  interference 
must  be  included  as  outlined  in  section  3.3-4.  Figures 
22-23  show  the  echo  train  as  photographed  on  the  oscillo¬ 
scope  screen  for  liquid  helium  at  4.1  K  and  for  seven 
crystals  near  melting. 

Two  difficulties  arise.  If  the  reflector  is 
tilted  then  the  beam  will  wander  further  in  the  -x 


* 


■ 


Figure  22  -  Echo  Trains  at  5  MHz 
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.80  K  (d)  Crystal  A6-e  3.18K 


Figure  23  -  Echo  Trains  at  5  MHz 
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direction  (as  in  figure  10(a))  with  each  reflection  until 
it  will  finally  disappear  for 

n  *  a/(2Dd)^  .  (4.2-1) 

At  5  MHz  this  will  occur  at 

n  *  7  n^  (4.2-2) 

max 

where  n  is  the  echo  number  of  the  first  maximum, 
max 

Typically  then  the  signal  should  disappear  for  n  v  40, 
but  it  obviously  does  not.  This  difficulty  is  avoided  if 
the  reflector  is  shaped  as  in  figure  10(b).  Secondly,  for 
a  tilted  reflector,  the  interval  from  n  =  0  to  the  first 
minimum  should  be  approximately  twice  as  long  as  from  the 
first  minimum  to  the  first  maximum.  Such  a  feature  is 
not  always  observed,  as  for  crystal  A6-b  in  figure  22. 
Again,  the  curved  reflector  is  less  restrictive  and  can 
fit  the  observed  envelope  reasonably  well  over  the  first 
two  maxima.  The  attenuation,  then,  is  analyzed  using 
equations  (3.3-32)  to  (3.3-42). 

-5 

3  in  equation  (3.3-24)  ranged  from  5x10  to 
10-4  cm-1,  corresponding  to  a  radius  of  curvature  of 
104  to  2xl04  cm.  This  compares  with  a  radius  of  ^  104  cm 
expected  from  the  elastic  modulus  of  beryllium  copper. 

It  is  interesting  to  note  that  the  attenuation  near  the 
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melting  point,  about  0.01  cm  ^ ,  is  only  4%  of  the  values 
of  around  0.25  cm  ^  that  one  would  obtain  by  looking 
only  at  the  echoes  previous  to  the  first  minimum,  as  was 
done  by  Hiki  and  Tsuruoka  (1976) . 

4 . 2-2  Attenuation  Results 

Figures  24-25  show  the  echo  train  at  low  tempera¬ 
tures  for  six  crystals,  as  well  as  the  15  MHz  and  25  MHz 
signals  for  crystals  A8-d  near  melting.  Attenuation  as  a 
function  of  temperature  and  frequency  is  illustrated  in 
figures  26-31. 

For  most  crystals  the  attenuation  follows  a  con¬ 
sistent  pattern.  At  higher  temperatures  it  rises  gradu¬ 
ally  with  decreasing  temperature  and  then  more  quickly 
as  the  onset  of  the  velocity  anomaly  is  approached.  It 
continues  to  rise  to  the  lowest  temperatures  measured. 
There  is  a  suggestion  of  an  attenuation  peak  in  some 
crystals  at  low  temperatures,  especially  at  15  MHz,  but 
the  large  measurement  uncertainties  preclude  any  definite 
conclusion.  The  attenuation  increases  by  as  much  as  100 
times  between  melting  and  0,3  K.  Most  interesting  are 
the  two  crystals  A6-b  and  A7-a.  The  absence  of  an  anomaly 
in  these  is  dramatically  confirmed  by  the  attenuation, 
which  increases  only  slightly  to  the  lowest  accessible 
temperatures.  At  5  MHz  most  crystals  which  show  100  or 
so  echoes  at  melting  have  only  2-4  remaining  at  0.3  K. 


Figure  24  -  Attenuation  of  Echo  Trains  at  5  MHz 
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Figure  25  -  Attenuation  at  5  MHz  and  Echo  Trains  at  15  MHz  (A8.e)  and  25  MHz  (A8.f) 
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Figure  26  Attenuation  as  a  function  of  temperature  and  frequency. 
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Figure  27  Attenuation  as  a  function  of  temperature  and  frequency. 
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Figure  28  Attenuation  as  a  function  of  temperature  and  frequency. 
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Figure  29  Attenuation  as  a  function  of  temperature  and  frequency. 
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Figure  30  Attenuation  as  a  function  of  temperature  and  frequency. 
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Figure  31  Attenuation  as  a  function  of  temperature  and  frequency. 
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These  two  exceptions  still  show  40-60  visible  echoes  at 
that  temperature. 

There  is  a  consistent  rise  in  attenuation  with 
increasing  frequency,  especially  between  5  and  15  MHz, 
where  a  is  as  much  as  20  times  larger  at  the  higher 
frequency.  The  25  MHz  results  generally  show  a  still 
higher  value,  although  the  uncertainties  are  rather  large 
at  this  frequency. 


CHAPTER  5 


DISCUSSION  OF  RESULTS 


5 . 1  Dispersion 

For  a  simple  linear  chain  model  one  would  only 

-9 

expect  a  velocity  change  of  Av/v  ^  10  between  5  and  15 
MHz.  The  observed  effect  is  five  orders  of  magnitude 
larger  and  so  cannot  continue  very  long  before  some 
anomalous  dispersion  sets  in  at  higher  frequencies.  In 
fact  two  crystals,  A8-a  and  A8-d,  do  already  show  such 
anomalous  dispersion  between  15  and  25  MHz.  We  cannot 
predict  the  extent  of  the  positive  dispersion,  particular¬ 
ly  if  it  ranges  all  the  way  up  to  thermal  energies  (in 

3 

this  experiment  kT/hu)  ^  10  )  .  If  this  were  so  then  fii<l 
to  much  lower  temperatures  than  expected  from  thermal 
measurements,  since  three  phonon  scattering  processes 
would  be  allowed  (section  2.6).  In  liquid  helium  such  an 
effect  has  been  described  by  Benin  ( 19 7 6 )  ,  Wehner  (1974) 
and  Maris  (1973) .  They  calculate  that  the  existence  of 
both  small-angle  three  phonon  scattering  and  wide-angle 
scattering  formed  from  many  back-to-back  small-angle 
collisions  result  in  a  region  of  second  sound  propagation 
existing  well  into  the  got>1  regime.  The  transition  to 
ballistic  propagation  also  becomes  more  complicated. 

There  is  evidence  of  wide-angle  scattering  in  solid  helium, 
as  discussed  by  Rogers  (1972) . 
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It  is  uncertain  from  the  data  if  the  large  dis¬ 
persion  is  in  some  way  linked  to  the  velocity  anomaly. 

On  the  one  hand,  the  crystals  which  exhibit  no  anomaly 
still  have  a  large  dispersion,  but  on  the  other  hand  it 
may  be  significant  that  the  two  crystals  which  show  the 
velocity  anomaly  beginning  to  disappear  at  low  tempera¬ 
tures  also  exhibit  positive  dispersion  at  25  MHz. 

5 . 2  The  Adiabatic  Region 

Equation  (2.5-10)  is  a  very  good  representation  of 
the  adiabatic  behaviour  of  the  sound  velocity.  Unfortun¬ 
ately  (2.5-9)  does  not  work  at  all.  Ahlers  (1970)  obtained 
the  Griineisen  constant  from  specific  heat  as: 

y  =  1.02  +  0.083  V.  (5.2-1) 

Then  [y  -  dtn  y/d£n  V]^  =  4  =  1-81  ,  (5.2-2) 

so  that  (2.6-10)  predicts  a  velocity  increasing  with 
rising  temperature.  Obviously  then  the  use  of  an  averaged 
Griineisen  constant  may  be  valid  for  thermal  measurements, 
but  it  breaks  down  completely  when  the  sound  velocity  is 
considered;  then  equation  (2.5-6)  must  be  used.  Unfor¬ 
tunately  almost  nothing  is  known  about  the  Griineisen 
constants  for  individual  modes. 


Because  the  constants  a  and  b  are  frequency 
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and  orientation  dependent,  the  assumption  of  a  reduced 
equation  of  state  which  led  to  (2.6-10)  is  invalid.  This 
conclusion  has  been  reached  previously  by  Jarvis  et  al 
(1968)  . 

5 . 3  The  Velocity  Anomaly 

From  figure  4  we  see  that  mx  ^1  for  5  MHz  sound 
in  the  temperature  range  1.5<T<1.8  K.  This  corresponds 
roughly  to  the  onset  temperatures  of  the  velocity  anomaly 
in  various  crystals.  However  figure  4  predicts  that  the 
anomaly  should  begin  at  a  10%  higher  temperature  for 
15  MHz  sound;  the  opposite  behaviour  is  observed.  The 
large  tn  does  not  permit  any  region  (dtn<1<ojtu,  although 
it  may  be  that  cotn<1  down  to  1.2  K  (Rogers  (1972),  Bert- 
man  et  al.  (1966),  Franck  and  Hewko  (1973)).  Still  there 
would  be  at  best  be  only  a  very  narrow  second  sound 
window  and  none  at  all  for  higher  frequencies,  unless 
there  is  some  interaction  with  the  "driftless"  type  of 
second  sound  which,  according  to  Hardy  (1970) ,  can  propa¬ 
gate  at  low  temperatures.  In  any  case  the  effect  predict- 

-4 

ed  by  equation  (2.6-10)  is  only  about  Av/v  ^10  at  1.2  K 

4 

and  decreases  as  T  . 

The  transition  to  zero  sound  is  more  complicated. 
In  equation  (2.7-3),  assume  an  isotropic  material  with  a 


unique  y  . 
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Then 


Y  _  2 

ef  f 

Y 


'K ♦ v  (kj) 

9 

n-K*v  (kj)/  p 
y 


If  the  dispersion  is  given  by 


(5.3-1) 


(jo  =  ck  [1-6  (k)  ]  ,  6  (k)  <<1, 


(5.3-2) 


then  (5.3-1)  becomes, 
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eff 


Y 


■  /lh±L 

~  \l-[l- 


—  6  (k) +6 (K) -k6 1  (k)+K6 1  (K) ] cos 0  \ 

6 (k) +6 (K) -kS ' (k) +K6 1 (K) IcosG/  P  * 


(5.3-3) 


Now  using  a  Debye  model  with  x  =  X  =  3"h^f 

and  after  performing  the  angular  integration. 
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e_f_f 

Y 


15 

811' 


-.  P 


[ink  |  6(x/Bftc)  -6  (X/3tic)  + 


+  x<5 '  (x/$fic)/$fic  -  X5 '  (X/3fic)/3fic  |  +  2] 


x  4  , 
e  x  dx 

,  x  ,  v  2 
(e  -1) 


(5.3-4) 


We  can  approximate  (5.3-4)  by  assuming  an  approximately 
constant  value  for  the  logarithmic  factor  (the  x  =  1 
value)  for  x  >  a""1  where  a  is  an  arbitrary  constant  such 
that: 


(8ftft)  1  >>  a  >  1. 


(5.3-5) 
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Then : 


eff 
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[ inh I  5 (x/3 he)  - 


'  o 

+  x6  1  (x/8T ic)/8hc  -  X6  '  (X/ghc) /O^c  | 

-  1  -  ^£n^[6(l/6ti c)  -6(X/8hc)  +  6' 

-  X6  '  (X/$fic)/6fic  J  . 


6 (X/3hc)  + 
+2]  x2  dx 
( 1/gfic)  /gfic 


(5.3-6) 


If  we  choose  simple  normal  dispersion. 


6 (k)  =  Ak2  , 


(5.3-7) 


then , 


From  (5. 
A  =  2x10 
observed 
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2  (3fic) 


(5.3-8) 


3-8),  (y  ff/y)2  %  2,5  for  T  =  1  K  and 

“16.  This  value  is  too  small  to  produce  the 


4 

anomaly,  even  above  1  K,  and  the  T  temperature 


dependence  would  make  it  negligible  below  0.5  K. 

A  reduction  in  the  normal  relaxation  time  could 
affect  the  zero  sound  regime  as  well  as  the  second  sound 
region,  but  again  the  anomaly  should  shift  to  a  higher 
temperature  for  a  higher  frequency. 

The  Granato-Liicke  dislocation  theory  (equation 


. 
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2.8-11)  can  be  made  to  fit  the  data  well  in  almost  all 

cases.  However  this  is  not  surprising  when  one  considers 

that  there  are  five  parameters  (A^,  A2,  k  and  g  and  n  in 

d  =  gTn)  that  can  be  varied.  In  fact  there  are  only  four 

independent  parameters  if  (2.8-13)  is  used.  But  for  an 

exponential  distribution,  k  >  2,  contrary  to  what  Wanner 

et  al.  (1976)  find  in  at  least  one  case. 

When  (2.8-11)  was  fitted  to  the  data,  n  varied 

-7 

between  2.7  and  5,  while  g  was  of  the  order  of  10 
2 

dyne-sec/cm  .  This  value  is  three  orders  of  magnitude 
larger  than  expected  from  Brailsford  (1971)  where  g  is 
given  as, 


g 


U5)3 

Trpb2 


R(0) 


0 


(5.3-9) 


where  £  is  the  Riemann  zeta  function,  8  is  a  numerical 
constant  -  1.2xl04,  b  is  the  Burgers  vector  and  R(0)  is 
given  in  terms  of  the  dislocation  contribution  to  the 
thermal  resistance  WD  as: 

R ( 0 )  =  WdT2/A|t=0  .  (5.3-10) 


Since  dislocations  generally  lie  in  the  basal  plane  of 
hexagonal  crystals,  especially  when  the  ratio  of  lattice 
parameters  c/a  is  the  ideal  spherical  atom  value,  (Hirth 
and  Lothe  (1968) ) ,  then  we  let  b  =  a. 


From  Klemens  (1958) , 


. 
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.2  2 

h  VDY 


423k' 


(5.3-11) 


The  g  ^  8x10  ^  dyne-sec/cm^ . 

A  was  of  the  order  of  10  cm  ,  an  order  of 

magnitude  smaller  than  Wanner  et  al,  (1976)  found.  For 

the  two  crystals  A6-b  and  A7-a  only  an  upper  bound  for  A 
2  -2 

of  10  cm  can  be  used.  We  can  estimate  the  dislocation 
density  from  the  mosaic  spread  (Hirsch  (1956))  seen  in 
neutron  scattering  experiments  (Minkiewicz  et  al»  (1973)). 
The  lower  bound  for  this  estimate  would  be  2x10  cm  , 
but  this  figure  is  a  gross  underestimation  because  it 
assumes  that  all  dislocations  bend  in  the  same  direction. 
An  upper  bound,  calculated  by  assuming  a  random  distribu¬ 
tion,  should  be  much  closer  to  the  actual  value;  this 

9  -2 

procedure  gives  A  ^  2x10  cm  .  I  conclude  that  A  is 
•  7  f!  “2 

m  the  range  10-10  cm  .  The  small  unphysical  A 
needed  for  a  fit  to  the  data  compensates  for  the 
unphysically  large  g  required. 

Although  almost  any  set  of  data  could  be  fitted 
well  on  its  own,  the  5  MHz  constants  did  not  give  a  good 
fit  at  15  MHz  or  25  MHz  for  the  same  sample.  The  magni¬ 
tude  of  the  observed  anomaly  did  not  decrease  as  fast  as 

-2 

the  approximately  oj  variation  predicted  by  the  theory. 
More  significant  was  the  temperature  of  the  anomaly, 
represented  by  the  onset  temperature  or  the  temperature 


’ 

X- 


Ill 


Tm  of  maximum  slope  d(Av)/dT.  From  (2.8-13),  TM  for  any 
one  crystal  is  given  quite  accurately  by 

TMn/w  =  constant.  (5.3-12) 

Equation  (5.3-12)  predicts  that  the  anomaly  should  occur 
at  a  temperature  20  to  50%  higher  at  15  MHz  than  at  5  MHz, 
but  in  all  cases  the  anomaly  shifts  to  a  temperature  15 
to  35%  lower.  This  argument  rules  as  decisively  against 
the  dislocation  theory  as  against  any  relaxation  time 
theory  discussed  previously  in  this  section.  In  fact  any 
resonance  type  of  theory  will  run  afoul  of  this  same 
objection  -  the  unusual  frequency  dependence  of  the 
anomaly . 

Finally  it  should  be  noted  that  although  (2.8-11) 
allows  for  a  negative-  or  positive-going  anomaly,  the 
negative-going  interval  must  occur  at  a  higher  temperature, 
so  that  once  d(Av)/dT<0,  it  cannot  be  positive  again  at  a 
lower  temperature.  This  predicted  behaviour  is  contra¬ 
dicted  by  crystal  A8-d. 

Saslow  (1977)  has  estimated  a  velocity  change  of 

-4 

Av/v-10  from  equation  (2.9-8).  Equation  (2.9-10)  can 
be  written: 

-r  T  /3V\  .1  /8T\  _  1  /3V _ \ 

0  1  V  \8T/  y,e  T  \dS J\ife  V  Kde^J  S,y 

(5.3-13) 
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(5.3-14) 


where  a  is  the  thermal  expansion  coefficient.  For  a 
perfect  crystal  the  bracketted  term  in  (5.3-14)  is  zero. 
If  Schottky  vacancies  with  energy  of  formation  (J)  are 
present  in  a  concentration  of 


n/N 


-<j>/kT 

9 


(5.3-15) 


from  Miyoshi  et  al.  (1970)  and  Reich  (1963)  , 


<J>  «  V 


(5.3-16) 


approximately  so  that 


d V 
V 


6d>  n  \  , 
kT* N  )  d£ll  * 


(5.3-17) 


Using  (5.3-17)  in  (5.3-14),  we  find  finally: 


t  aT  6_^_  -<J>/kT 

0  ~  3  kT  e 


(5.3-18) 


-5  3 

From  Jarivs  et  al. (1968) ,  a  is  approximately  6x10  T 
K-*1  at  17.4  cm3/mole.  <j>  is  50  ±10K,  so  that: 


6 


2x10 


300 

T 


e-SQ/T 


(5.3-19) 
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At  2K,  6  =  3x10  4  and  at  1  K  is  2x10  5.  The  large  volume 
dependence  of  <j>  more  or  less  rules  out  vacancies  as 
making  a  significant  contribution  to  the  crystal  proper¬ 
ties  at  this  density.  As  with  zero  sound,  the  T4  temper¬ 
ature  dependence  of  the  effect  should  obliterate  it  below 
1  K.  Also  c^  /ci  <  1/  so  that  the  effect  is  really 
rather  small  (and  negative) ,  although  the  fact  that  c^2 
rises  with  decreasing  temperature  could  moderate  the 

4 

T  dependence  in  (5.3-19). 

In  general  the  velocity  anomaly  and  other  observed 
effects  follow  no  discernible  pattern  from  crystal  to 
crystal.  There  does  not  appear  to  be  any  dependence  upon 
orientation.  Although  the  crystals  were  grown  all  in  the 
same  manner,  one  is  led  to  examine  defect  mechanisms  since 
defect  types  and  concentrations  can  vary.  The  range  of 
dislocation  densities  found  from  the  Granato-Luche  theory 
appears  too  low.  In  equilibrium  the  vacancy  concentration 
is  much  too  low.  One  cannot  however  rule  out  other  non¬ 
equilibrium  defects  such  as  clusters  of  vacancies  or  stack¬ 
ing  faults.  Unfortunately  almost  nothing  is  known  about 
these  structures  in  solid  helium. 


. 
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5 . 4  Attenuation 

From  equation  (2.5-5),  the  attenuation  in  the 
adiabatic  region  will  be 

a  <  4xl0~8  T  e16 • 6/T  cm_1,  (5.4-1) 

or  2x10  ^  cm  ^  at  3.5  K  and  7x10  ^  cm  ^  at  1.8  K.  It 
would  of  course  be  impossible  to  measure  such  small 
attenuations.  The  difficulty  of  matching  impedances  to 
obtain  a  good  signal  for  the  higher  harmonics  must  almost 
certainly  have  contributed  to  the  measured  attenuation 
increasing  faster  than  linearly  with  frequency. 

In  a  hypothetical  second  sound  region,  the  attenu¬ 
ation  from  equation  (2.6-9)  and  figure  4  would  be: 

a  -  8x10  ^  T^/t  .  (5.4-2) 

At  T  =  1  K,  a  then  would  be  0.25  cm  ^  or  more,  which  is 
near  the  measured  values.  However  this  attenuation  should 
drop  with  decreasing  temperature  as  far  as  the  Poiseuille 
flow  maximum.  The  attenuation  is  not  frequency  dependent. 

Generally  all  theories  fail  because  of  the 
continuation  of  the  high  attenuation  over  such  a  large 
temperature  range.  For  the  dislocation  theory  also 
(from  (2.8-1)),  when  d<oo, 
a  oc  d  a  Tn  , 


(5.4-3) 


. 
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so  that  the  attenuation  should  be  falling  rapidly  by  the 
time  the  plateau  in  the  velocity  is  reached.  This  situa¬ 
tion  does  not  occur.  One  expects  that  there  will  eventu¬ 
ally  be  an  attenuation  peak  and  a  0  as  T  ->  0 ,  but  this 
behaviour  has  not  yet  been  observed. 


CHAPTER  6 


CONCLUSIONS 


There  is  a  very  large  frequency  dispersion  at  all 

,  4 

temperatures  in  pure  single  crystals  of  hep  solid  He  . 

The  dispersion  is  usually  normal,  but  both  direct  and 

indirect  evidence  exists  for  anomalous  dispersion  at  or 

above  25  MHz.  The  longitudinal  sound  velocity  rises  with 

decreasing  temperature  until  an  anomaly  sets  in  between 

1  and  1.8  K.  The  velocity  changes  very  rapidly  with 

temperature  over  a  range  of  0.5  -  1  K  and  then  levels  off 

_3 

to  a  constant  value  with  (v/v  1  ^  1-3  xlO 

aaiaDatic 

Two  crystals  show  a  deviation  from  this  plateau  below 
about  0.3  K.  In  most  samples  the  anomaly  is  positive, 
but  two  crystals  show  a  negative  anomaly  while  two  others 
exhibit  no  anomalous  behaviour  at  all.  The  anomaly 
becomes  smaller  and  moves  to  lower  temperatures  as  the 
frequency  is  increased. 

The  attenuation  is  very  small  in  solid  helium 
near  melting,  generally  0.01  cm  1  or  less.  Very  careful 
consideration  must  be  given  to  the  echo  envelope  when 
measuring  these  low  attenuations.  The  attenuation  rises 
as  the  temperature  is  reduced,  slowly  at  first  and  then 
more  quickly  as  the  velocity  anomaly  is  reached.  The 
attenuation  continues  to  rise  to  at  least  0.25  K  or  lower 
There  is  no  theory  at  present  which  can  explain 
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the  behaviour  of  the  velocity  and  attenuation  in  the 

anomalous  region.  All  theories  fail  completely  on  two 

counts:  the  anomaly  moves  to  lower  temperatures  as  the 

frequency  is  increased;  the  attenuation  continues  to  be 

very  large  down  to  0.25  K. 

Future  work  should  proceed  in  three  directions : 

to  lower  frequencies,  to  lower  temperatures  and  to  an 

investigation  of  transverse  sound  propagation.  Most 

useful  of  all  would  be  a  similar  experiment  done  on  very 
3 

pure  He  ,  in  both  the  hep  phase  (to  see  if  isotopic 
effects  are  important)  and  the  bcc  phase  (to  observe  the 
significance  of  the  crystal  structure) . 
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APPENDICES 


A1  Extrapolation  of  Elastic  Constants 

Data  for  this  extrapolation  were  taken  from  Wanner 
and  Franck  (1970),  Greywall  (1971)  and  Crepeau  et  al. 

(1971)  . 

There  are  several  useful  relations  among  the 
elastic  constants.  From  Franck  and  Wanner  (1970), 


C11  +  c12  c13  +  c33‘ 


(Al-1) 


If  (Al-1)  is  used  in  the  expression  for  compressibility 
(Nye  (1957)),  then 


kt  =  3/  (C33  +  2c13)  .  (Al-2) 

c^  was  fitted  to  an  empirical  relation  using  all  data 
(but  especially  those  of  Greywall  which  were  spread  over 
a  wide  volume  range)  to  give: 

c±1  =  4 . 236x10  X1e  ”v/3*017  .  (Al-3) 

Then  for  almost  all  data,  c^  and  c^  were  found  on 
average  to  be  simply  related  according  to: 


Cll//c12  =  ±  0*02. 


(Al-4 ) 


125 
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From  equations  (Al-1)  to  (Al-4) ,  the  four  elastic  constants 

cll'  C12'  c13  and  c33  can  be  calculated  if  kt  is  known. 

For  this  purpose  the  data  of  Jarvis  et  al.  (1968)  was  found 
to  be  fitted  very  well  by 


4.2513x10  12 


eV/3.0703 


(Al-5 ) 


In  order  to  calculate  c^ ,  we  look  at  the  Debye  velocity 
vD  given  by 


[  /  (v  3  +  V  3  +  v  3)  dn/12TT]'1/3, 

1  2 


as  in  equation  (2.4-8).  Now, 


(Al-6) 


(Al-7) 


and  @o  is  given  empirically  by  Ahlers  (1970)  as: 

0o(V)  =  68.56  (V/14.208)'1-02  exp  [-0 . 083 (V-14 . 208) ] . 

(Al-8) 

Then ,  using  equations  (2.3-9)  to  (2.3-12)  for  the  sound 

velocities,  c^4  can  be  calculated  on  a  computer  once  the 

other  four  elastic  constants  are  known. 

I  could  place  a  great  deal  of  confidence  in  the 

3 

extrapolation  to  17.4  cm  /mole  because  extending  the 

3 

extrapolation  all  the  way  to  16  cm  /mole  gave  a  value  for 
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c33  only  7%  higher  than  that  found  by  Reese  et  al  (1973) 
from  neutron  scattering  (25.4xl08  dyne/cm2). 

Some  representative  values  for  the  elastic 
constants  are  given  in  table  2. 

Extrapolation  of  Relaxation  Times 

From  the  measurements  of  Hogan  et  al.  (1969)  at 
85  and  126  atmospheres  and  the  more  accurate  measurements 
of  Lawson  and  Fairbank  (1973)  at  85  atm.,  the  normal 
relaxation  time  is 


TN  (2.02±0.1)  x  10  ^2  (0^/T) 2  sec. 


(A2-1 ) 

The  thermal  conductivity  in  the  umklapp  scattering  region 
was  fitted  accurately  by  the  expressions: 


K 


x 


K 


V 


Vb//  T  ' 


(A2-2 ) 


with  bx  =  2.58  and  b ^  =  4.62.  The  coefficients  Ax 

and  A  tl  are  functions  of  pressure.  The  data  is  scant  in 

this  respect.  However  Hogan  et  al.do  find  that 

Ax  (85  atm)  =  112  erg/cm-sec-K,  A^(85)  =  2310  erg/cm-sec-K, 

A„  (53.4)  =  0.61  A  a  (85)  and  Ax(126)  *  1.5  Ax(85). 

We  can  estimate  that  for  a  melting  pressure  P  and 
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TABLE  2 

ELASTIC  CONSTANTS 


Volume 

3 

(cm  /mole) 

cn 

C12 

(all  in 

C13 

units  of 

C33 

O 

10  dyne/cm 

c44 

2> 

21 

4.01 

2.25 

1.29 

4.98 

1.59 

20.5 

4.74 

2 . 66 

1.49 

5.91 

1.74 

20 

5.59 

3.14 

1.73 

7.00 

1.92 

19.5 

6.60 

3.71 

2.01 

8.30 

2.11 

19 

7.79 

4.38 

2.32 

9.84 

2.34 

18.5 

9.19 

5.17 

2.69 

11.66 

2.61 

18 

10.85 

6.10 

3.12 

13.83 

2.96 

17.5 

12.81 

7.20 

3.62 

16.39 

3.33 

17.4 

13.24 

7.44 

3.72 

16.95 

3.43 

17 

15.12 

8.49 

4.19 

19.42 

3.77 

16 

21.06 

11 . 83 

5.61 

27.28 

4.93 
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propagation  either  parallel  or  perpendicular  to  the 
c-axis , 


A (P)  =  (P/85)1,05  A (85) 


(A2-3) 


Now, 


x  =  3kV/Cvvd"  , 


(A2-4 ) 


so  we  have  finally: 


U 


1.73x10 


-8 


D 


V 


0 


D 

T 


ea 


1.05  0d/2.58T 

e 


(A2-5) 


U 


3.57x10 


-7 


v 


D 


V 


0 


D 

T 


(a 


1.05  0d/4.62T 

e 


(A2-6 ) 


Equations  (A2-1) ,  A2-5)  and  (A2-6)  are  plotted  on 
figure  4. 

The  total  relaxation  time  t  is  calculated  from 
the  following  relations  of  Hogan  et  al.  and  Guyer  and 
Krumhansl  (1966): 

(a)  Boundary  scattering  region  (vDxN>>d,  vDxu>>d) 


x  =  ad/v^  where  a  >  1 . 


(A2-7) 


2  2 

(b)  Poiseuille  flow  region  (vQxN<<d,  xu>>d  /v^  xN) 
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T  =  §2  d2/vD2TN-  (A2-8) 

(c)  Umklapp  region  (t  <<t  ,  xu<<d) 

T  =  Ta  .  (A2-9 ) 

The  total  relaxation  time  t  is  shown  for  a  cell  width 
d  =  1  cm  by  the  dashed  lines  in  figure  4  after  interpol¬ 
ating  equations (A2-7)  -  (A2-9). 


